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SUMMARY 


We describe a procedure for the generation and 
reduction of the radiometric data Known as 'REGRES' files. 
These data are a current output of NASA's Orbit 
Determination Program. The software package was developped 
in view of the data analysis of the gravitational wave 
experiment, (GWE) planned for the European spacecraft 
'Ulysses'. 



INTRODUCTION 


Doppler tracking of interplanetary spacecraft is 
routinely performed by the stations of NASA's Deep Space 
Network (DSN). An electromagnetic carrier of highly stable 
frequency is sent from a ground station to a 
spacecraft, which retransmits it back coherently to the 
Earth by means of a trasponder. 

The presently used radio link configurations for 
precision Doppler tracking are of two types: in the first 
one, called 'two-way' operation, the same station both 
transmits and receives the radio signal, in the second 
one, called 'three-way', the transmitting and receiving 
station are different. The signal is usually transmitted 
in S band (2.1 GHz) and received both in S and X band (8.4 
GHz). 

We call Doppler residual the quantity obtained by 
subtracting the measured frequency shift of the carrier to 
the predicted one. In the REGRES data, the predictions 
are obtained a posteriori by means of a regression 
method: the received radiometric data are used to improve 
the spacecraft ephemeris in order to minimise the least 
square deviation between the observables and the 
predictions. 

In this paper we report a description of the REGRES 
data generation and of the software we have developped to 
reduce and display the data. This software has been 
applied to the analysis of four days (DOY 308 - DOY 313, 
1980) of REGRES data relative to the Voyager I spacecraft. 

The software package was implemented on an HP-1000F 
computer, with an RTE-6-VM operative system. The graphical 
unit was an HP-9872T plotter. A general purpose graphical 
routine, developed at IFSI-CNR, was used to plot output 
data. 
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REGRES FILE GENERATION 


REGRES files are a current output of the Orbit 
Determination Program (ODP). The volume and the structure 
of REGRES files are not fixed, but depend on the way 
the ODP is run. Standard references on the subject are 
(Moyer, 1971) and (Khatib et al., 1972). ODP is actually 

a set of many links, which, to some extent, may be 

run independently. Its functions are integrated by a 
number of allied programs, the most important of which 
is DPTRAJ. The process which leads to the computation 

of the trajectory of the spacecraft is described in the 
following. 

The fundamental quantities handled by the ODP are: 

X=(r,r )'• spacecraft state vector, i.e. position and 

velocity of the spacecraft with respect to a center 
of integration (usually the baricenter of a celestial 

body). 

g = (X 0 ,a) : parameter vector, containing the initial state of 

the spacecraft (X 0 ) and a set (a) of physical parameters 
(see below) needed for the integration of the equations of 
motion. 

z: vector of observable quantities (Doppler, range, 

angle, ADOR). 

The main physical parameters a are: 

- the gravitational parameters p^GM* of Sun, 
planets and satellites of the solar system; 

- the coefficients of the harmonic expansion of the 
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planets' gravitational potential; 

- the quantities used in modelling the probe 
acceleration due to small forces and manoeuvres; 

- atomic station time -> UTC conversion parameter; 

- coordinates of the tracking stations; 

- the relativistic parameter y. 

Many other quantities may be used, depending on the 
particular spacecraft: there are parameters relating to 

atmospheric models, mass anomalies, gas leaks, solar 

radiation pressure, etc. 

The vector a is usually stored as a part of a file 
called Generalised Input File (GIN file), characterising 
the physical world in which the probe moves. 

The ODP not only performs the integration of the 
equations of motion, but also allows the 

redeterminations of the parameters g and, therefore, the 
determination of a better state vector X(t). This goal is 
obtained by comparing predictions and observations in a 
regressive analysis made up of several steps. 

STEP 1 PATH- VARY (or PVDRIVE) 


This link has actually a twofold purpose and may be 
separated for clarity into two sub-units (PATH and VARY). 
PATH integrates the equations of motion of the probe, 
starting from given injection parameters Xo and physical 
parameters a. The probe acceleration results from the 
following contributions: 

- Newtonian gravitational forces (leading term) from 
the relevant celestial bodies. 

- Oblatenesses of planets. 

- Mascons. 

- Relativistic effects. 

- Solar radiation pressure. 

- Attitude control. 
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- Orbital manoeuvers. 

- Atmospheric drag (in case of planetary encounter). 

The center of integration may he any planet, or 
center of mass of planetary system (i.e. planet with 
satellites). The output of PATH is a propagation of the 

initial state vector X along the trajectory as a function 
of time and of the parameters g: 

X : X ( t, 3) 


Its output X(t)=(r(t),r(t)) may be stored in the Probe 
Ephemeris Tape (PET). It is also used as input to the 
program PREDIX which provides observable predictions to the 
tracking stations. These predictions are used to generate 
the pseudoresiduals contained in the Archival Tracking Data 
Files (ATDF). 

The other subunit (VARY) integrates the variational 
equations 
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whose solution gives Z and Z as functions of time. 
These quantities are needed in the subsequent step 
(REGRES) to get the partial derivatives of the observable 
quantities with respect to the parameters. The matrices A, 
B and C are obtained from the equations of celestial 

mechanics. 


7 



STEP 2 (REGRES) 


In this step, the regression partial derivatives of 
the observables (Doppler, range, angle) with respect to 
the parameters g are formed. The regression partial 
derivatives are needed to adjust the quantities xo and a 
in order to obtain a better estimate of the spacecraft 
trajectory (i.e. better residuals). 

To this end, for each acquired value of the 
observable r at the reception time t 3 , the corresponding 
ground transmission time t t and spacecraft reception time 
t 3 are computed. The times tj and tg are obtained by 
solving the light time problem in the metric of the solar 
system baricenter. For each leg (uplink and downlink) the 
transmission and reception times are given by (neglecting 
terms of order i/c 5 ): 


t — + ^ 


r 1 J * + Y 

ps In 

(ritrj+rij) 

v re c urans - 

c 

(ri+rj-rij) 

where 





r i 1 i 

I Ii 

‘Vans’ II 



r j * i 

i 

‘Vec’ II 



V I 

i ^ 

(t ) - r 

rec _i 

(t 

trans 

> II 


P s r GM S = gravitational parameter of the Sun 

y = post newtonian parameter (in general 
relativity y = i) 


The vectors r have origins in the solar system 
baricenter. In the last version of REGRES bending terms as 
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well as Jupiter and Saturn terms are included. 

The solution of the light time problem must be 
performed for each acquired point, and requires therefore 
considerable computer time. 

The regression derivatives dr/dg are then 
formed by considering that the observable r is a 
function of the state vectors of the Earth stations 
and the spacecraft; 

Z = Z[X st (t 3 ,3) > X sc (t 2 ,3),X st (t 1 ,g);3) 

It may be seen that the partial derivatives 
dr/dg are actually the sum of many terms, originating 
from explicit and implicit dependence of Z on the 
parameter vector. We refer to (Moyer, 1982) for a 

detailed discussion of the subject. Here we just point 
out that the solution of the variational equations, 
which provides the quantity 8X sc /dg, is required for the 
computation of dZ/dg. 

A typical REGRES output file contains: 

a) values of the parameters g for which partlals have 
been computed 

b) the selected observables and their residuals 

c) the selected partial derivatives dr/dg. 

It must be noted that at this step the parameters 
g are divided into two groups: "solve-for" parameters, for 
which new values are computed; "consider" parameters, 
which are not corrected but whose errors are taken 
into account in the computation of the errors for 
the "solve-for" parameters. 


STEP 3 (ACCUME-SOLVE) 

These two links use the partials obtained from the 
previous step to compute the parameter vector g and 
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its covariance matrix. With the new q's, a new GIN 
file may be created, and a new iteration, PV, REGRES, 
ACCUME, SOLVE (PVRAS), may be started. 

The REGRES file which we analysed has not been 
generated by the complete sequence PVRAS, which is 
complicated and unnecessary for our purposes. The link 
REGRES was run in the so called "simulation mode", 
i.e. using already existing PV and GIN files and 
bypassing the computation of the partial derivatives ez/dg 
The following steps have been performed to get the final 

REGRES file: 

a) select the GIN file to be used, which contains an 
initial state vector (ro.JTo) of the spacecraft and a set 

of physical parameters a that already gave satisfactory 
residuals. 

b) select the corresponding PV file (containing 
probe ephemeris) 

c) select the file (ATDF or IDR) containing the 

observable quantities for which the residuals have to 
be computed (we actually used the same ATDF which was 

previously analysed at Frascati). 

d) run the program STRIPPER which, starting from an 

ATDF or IDR, generates an OD file (a file which can be 

used as an input to REGRES) 

e) run REGRES in the "simulation mode", using PV, GIN 
and OD files as input (together with the ephemeris of the 
bodies of the solar system). 

The output of the step (e) is the required REGRES 

file. 
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SOFTWARE 


In this section we report on the software package 
implemented, on the HPiOOOF System at CNR-Frascati to handle 
the REGRES files. 

These data, generated by an Univac 1108 System at 
JPL and stored on magnetic tape, have been converted 
using the Univac computer of the University of Rome in 

ASCII characters files stored on magnetic tape. The list 
of this Fortran program (named UNIHP) is shown in 
appendix B. 

The main program implemented on the HPIOOOF 
(named DECODEREGRESS) decodes the REGRES data (ASCII 
characters) and selects the relevant items for the GWE 
data. The format of the output data (called STANDARD format 
throughout the paper) is described in (less and Armstrong 
1985). 

The software package is completed by a set of 
programs for the display and validation of the data and 
for the storage of Doppler residuals (our observable) in 
the HP disk memory. These programs are named: 

STATREGRESS, PLOTREGRESS and FFREGRESS. 

STATREGRESS reads the STANDARD tapes and produces an 
output list for each day of data. This list contains 
hourly and total number of samples for each active 
station, band and mode. 

PLOTREGRESS provides plots of the STANDARD tapes. 

In order to simplify the access to the data, the 
STANDARD files can be copied on a disk using the FFREGRES 

program. Each disk file contains the data relative to a 

single pass. 

The program listings are reported in Appendix B. 
In the following we report a description of the software. 


11 



DECODEREGRES 


DECODEREGRESS reads the REGRES tapes In ASCII 
format, selects the items relevant to the GWE and generates 
the STANDARD tapes. 

DATA STRUCTURE 

INPUT 

BUFFER DIMENSION: ISK(40, 100) 

RECORD LENGTH: 6000 characters 

RECORD FORMAT: Each data record is composed of 100 cards 

of 80 characters each. The 80 characters of each card have 
the following meaning: 


char. 

1-5 : 

measurement group identifier Nl, ranging 



from to 99999 

char. 

6-8 : 

card identifier N2, ranging from l to 100 

char. 

9-32 : 

data value, written in the format F24. 18 

char. 

33-56: 

as above 

char. 

57-80: 

as above 


DATA FORMAT: In the original tapes generated at JPL, 

twenty quantities were associated to each Doppler data 
acquisition (see tab. 1). We will refer to this set of 20 
quantities as to a "measurement". Therefore each record 

contains 15 measurements (i. e.: 3*100/20), (see tab. 

2). A data block is defined as a set of 40 measurements. 
The block identifier number Nl, which appear in tab. 2, 
refers to the record generated at JPL with the Univac 
routine NTRAN. 

See tab. 3 for an example of ISK (40,100) record. 

OUTPUT 

The output data, stored on magnetic tapes, are 

structured as a sequence of fixed lenglft. records with eof. 
BUFFER DIMENSION: LMAT(160,15) 
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RECORD LENGTH: 2400 double integers number, i.e. 4800 

16-bit words. 

RECORD FORMAT: Each of the 160 rows contains 15 double 

integer items. 

DATA FORMAT: Each row contains information about one 
measurement in both bands (S and X). Tab. 4 shows the data 

structure of the STANDARD records. 

In tab. 5 we report, the items of the REGRES 
tapes generated at JPL and their corresponding items in 
the STANDARD records. Tab. 6 shows an example of 
STANDARD record (they e.c. refers to the same REGRES data 
in tab. 3). 

ALGORITHM 

The program reads in sequence the ISK(40,100) 
records. The 20 quantities corresponding to a single 
measurement are transferred in the INTEGER*4 buffer 
ISK1(240). This buffer is decoded and the relevant items of 
each measurement are transferred in a second INTEGERS 
buffer LMAT(160,15). It is important to note that two 
subsequent measurements in input usually refers to the S 
and X band respectively. In order to properly fill the 
eighth column of the matrix LMAT (S-3/11X), the program 
requires that S-band datum must precede the X-band datum in 
the input data. In fig. 2a and 2b is reported the 
flowchart of the program. 

The program is structured as a main routine, 
managing the input data buffer, and a subroutine, managing 
the output buffer. Those two program segments are now 
described in more detail, referring to fig. 2a and 2b. 


MAIN PROGRAM 

1) INFORMATION FOR INPUT/OUTPUT 

la) options and parameters for the Job: 

sequential number of the required initial REGRES 
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record, sequential number of the required final 

REGRES record; 

tape logical unit number; 

etc. 

lb) search for initial record to analyse 

2) READ A ISK (40. 100) RECORD 

this step is made using the system routine XTAPE 
(EXEC) 

check on the record length 
eof causes jump to step #8 
read error causes jump to next record 

3) TEST FOR COLUMN IDENTIFIER 

if the block identifier is changed, the program 
resets the decodifying buffer pointer IP. This step 
is needed to eliminate the free item present at 
every changing group 

4) DECODIFICATION AND STORAGE OF THE DATA 

this step is performed by the subroutine CONTROLLO 
which transfers and decodes the data 

5, 6, 7) REPEAT STEPS 1 TO 4 FOR ALL REQUIRED RECORDS 
AND TAPES 

8) FINAL OPERATIONS 

8a) stores the last data in memory 

8b) writes eof on STANDARD tape and rewinds the tapes 
8c) prints information on the work just carried out 

9) STOP 
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STATREGRESS 


THIS program reads the STANDARD tapes and produces an 
output list for each day. The list contains hourly and 
total number of samples for each active station, hand and 
mode. An example of an output list is reported in tab. 7. 

ALGORITHM 

The input data are contained, in the STANDARD files 
generated by DECODEREGRESS. The algorithm uses a buffer of 
sample counters structured as a matrix 48*24. Each row 
represents one of the combinations of 6 DSN, stations 
4 tracking modes and 2 bands, while the 24 columns 
represent the hours of the day. The identification of the 
sample counter for a given DSN station, band, mode and hour 
of day combination is given by: 

ROW = ( S — 1 ) *8+ (M-i) *2+B 

COL - HOUR OF DAY 

where B is i or 2 for S or X band respectively M 
represent the 4 different radio link modes, while 
S= 1, 2, 3, 4, 5, 6 represents the six DSS: 12, 14, 42, 43, 61 

and 63 respectively. The decodification of the row number 
can be made as follows: 

STATION = INT ( ( ROW - 1 ) / 8 ) + 1 

BAND = ROW'S PARITY 

MODE = INT ( ( ROW - BAND - ( STATION - 1 ) * 8 ) / 2 ) + 1 

In fig. 3 is reported the flow chart of the program, whose 
main steps are described below. 


1) ASK INFORMATION FOR INPUT/OUTPUT 
tape logical unit number 
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number of records to analyse 
Initial and final day and year 
etc. 

2) READ A RECORD 14 (160, 15) 

3, 4) SEARCH FOR FIRST RELEVANT DATA 


5,6,7) ANALYSE DATA AND FILL THE LIST 
6) OUTPUT OPERATIONS 

if the list (i. e. the output buffer) of the required 
day is empty, (if there are no data, for example) 
the program prints a warning, otherwise it prints 
the required list. 

10) START ANALYSIS OF A NEW REQUIRED DAY (if any) 

11) STOP 
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PLOTREGRESS 


Tills program provides plots of STANDARD data for eacli day 
and mode. 

DATA STRUCTURE 

The input data are the STANDARD file generated by the 
program DECODEREGRESS. The user selects the passes to be 
plotted for each desired day. DQY, mode and the DSS are 
required as input parameters. The data appear into 
three curves (S-band, X-band and S-3/il#X) as shown in fig. 
4. 


ALGORITHM 

The program flow chart is shown in fig. 5. 

1) ASK INFORMATION FOR INPUT/OUTPUT 

ia) options and parameters for the job as: 
initial and final REGRES records 
tape logical unit number 

required days and modes 
etc. 

ib) search for the initial record to analyse 

2) REWIND TAPE 

3) READ A RECORD 

4) CHECK FOR DAY AND MODE 

if data refers to the required day and mode, they 
are stored in a buffer 

5) DO DATA REFER TO A LATER DAY? 

in this case the selected data are plotted 

6) LAST RECORD? 

the search for the relevant data is repeated for 
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all records; 

the search ends if the data of a record refer to a 
later day; 

If there are no data to plot, the program prints a 
warning and goes to step #8. 

7) PLOTTING 

this step is carried out by a software package 
deve lopped at IFSI 

8) OTHER PLOTS? 

in this case, steps (1-7) are repeated. 

9) STOP 
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FFREGRES 


This program stores S and. X hand Doppler residuals 
and oscillator frequency from STANDARD tapes into disk 
files. Each file contains data referring to a given pass. 
Besides Doppler residuals, each file contains a header 
which reports spacecraft angular data, RTLT and other 
information. Some of these quantities are retrieved from 
the file DATA/FILE#i (see less and Armstrong, 1985) 


DATA STRUCTURE 

INPUT : The input data are the STANDARD data and 
angular data FILE#i (see Tab. 8) 

OUTPUT: In output one obtains ASCII disk files 

relative to a single pass 
BUFFER DIMENSION: I STORE ( 3, 730) 

RECORI^ENGTH : 2190 double integers, i. e. 4380 
16-bit word 

RECORD FORMAT: each of the 3 rows contains 730 
double integers. The file name has the 

following mnemonic code: 

DAT/RODTTRRDDD 

where : 

DAT is the data subdirectory; 

ROD stands for REGRES Original Doppler resid. ; 

TT is the code for transmitting DSS; 

RR is the code for receiving DSS; 

DDD are three characters indicating the DOY. 

DATA FORMAT: each file (see tab. 9) is composed by a 

header, with general information about the pass, 
followed by S and X band residuals and oscillator 
frequency in the first, second and third row 
respectively. Residuals are given in mmHz while 



the oscillator frequency is given in Hz. The header 
contains initial and final time, number of data, 
receiving and transmitting DSS, tracking mode, 
RTLT, spacecraft's right ascension and declination 
at the middle of the pass. An example of residual 
file is reported in tab. 10. 

ALGORITHM 

The algorithm looks for data referring to the 
required receiving DSS, mode and begin of tracking day. 
Then it selects data from FILE#l to build the header (see 
fig. 6). 

For the Voyager I data the sampling rate was of 1 
pt/min: therefore each Doppler residual refers to a given 
minute of the day (MOD). The time T (in MOD) of every 
datum can be obtained from the following relation: 

T = module ( I - li - T in , 1440) 

where I is the column number of the datum in the array 
ISTORE, and T^ n the MOD of the first datum of the file. 
Since the header must be skipped, 11 must be 
subtracted from I to get the right position. Finally, 
the remainder of the division of I - 11 -Tin by 1440 
provides the MOD of the datum (if I-ii-Tin > 1440, 
the MOD refers to the next day). 
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APPENDIX A 


LIST OF ACRONYMS AND NAMES 


ATDF 

Archival Tracking Data File 

DOY 

Day Of the Year 

DSN 

Deep Space Network 

DSS 

Deep Space Station 

EOF 

End Of File 

GIN 

Generalised INput file 

GWE 

Gravitational Wave Experiment 

I DR 

Intermediate Data Records 

HOD 

Minute Of the Day 

ODP 

Orhit Determination Program 

PET 

Prohe Ephemeris Tape 

PV 

Path-Vary 

REGRES 

Regression files 

RTLT 

Round Trip Light Time 
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APPENDIX B 


PROGRAM LISTINGS 


We report In this appendix the fortran program 
listings of the computer programs used in the data 
analysis. 
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nnn nnnnnnnn 


PROGRAM UN I HP 



tt*****-*******^™^^^^^ 7 ^ -' '• . ,„,_ 

* TRA5F0RMA7 1 ONE DEI NRbTRT REGRESS IN FORMRTO UNI UPC 

***^***X*v4***********************^ " — 

. 

* QUESTO PROGRRMMfl GPERP SUI DPT I REGREbb 

* 

**.***********-* 


* 

-* 

-X 


DICHIRRRZIUNI _ T 

COMMON / DPT I / R C SOI ) , ROUT (100), NROUT 
DOUBLE PRECISION PRRPY ( 20 , 40 ) , R 
CHRRRCTER*8G ROUT 


I0UT=10 
I0UT2 = ri 
PR I NT *1 

1 FORMPT( 2X , ' N . RECORD INITIAL AND FINRL?') 
RERDC5,*) NI N , NFI 

NRIN=NFI -NIN+1 
PR I NT 2 

2 F0RMRTC2X , ' DOYOU WANT PRINTOUT OF RC20,40)? 
READ (5,*) ISTRM 

NR = 0 
NROUT = 0 
LL = 0 
N 1 = N I N - 1 

IFCN'I . EQ. 0) GO TO 100 
200 CONTINUE 

CRLL NTRRNS ( I OUT, 2 , 1600, RRRRY , L ,22) 

nr=nr+-i 

IFCL.LT. 0) GO TO 500 
IFCNR.LT. NI ) GO TO 200 
100 CONTINUE 

CRLL NTRRNS C I OUT ,2, 1600 , RRRRY , L , 22 ) 

nr=nr+i 

IFCL.LT. 0) 60 TO 500 
IFCL.EQ.-2) GO TO 610 
IFC ISTRM. NE. 1 ) 60 TO 150 
WRITE C 8 , 300) NR 
WRITE (8 , 400) 

WRITE C 8, 700) (RRRRY CL , 1 ) ,L=1 , 20) 

WRITEC8, 400) 

WRITE (8,700) C RRRRY CL, 40) ,L = '’ ,20) 

WRITEC8 , 400) 

150 CONTINUE 

RC801 ) =0 . DO 
11=0 

DO 10 1=1,40 
DO 10 J=1 ,20 
11=11+1 

RC I I ) = RRRRY ( J, I) 


<1=51 >' ) 



u u 


10 CONTINUE 

DO 20 L=1 ,801 ,3 
LI =L+2 

COLL FILBUFCLL, I0UT2) 

ENCODE C 80 , 800 , AOUT (LL ) , JJ , ERR = 900) NR , LL , COCK) , K=L , LI ) 

IFC I5TAM . NE . 1 ) GO TO 20 

IFCL'I .GT. 9. AND. LI .LT.793) GO TO 20 

WRITE C 8 , 850 ) LL , ROUT ( LL ) 

850 FORMAT C 2X , ’ LL= ' , I3,2X,A80) 

20 CONTINUE 

600 FORMAT (I 5, 13,3024.18) 

350 FORMAT C2X, 'ERROR WRITE K , LL , NR , NROUT ' , 413) 

C 

IFCNR.GT.NFI ) GO TO GOO 
GO TO 100 
C 

300 FORMAT C 2X , ' RECORD NUMBER =',110) 

400 FORMAT ( / / ) 

700 FORMAT ( 2X , 3D24 . 18) 

C 

500 CONTINUE 

PRINT 3, NR, NROUT 

3 FORM AT C2X, 'HARDWARE ERROR, NR, L NROUT' , 2X,3I9) 

STOP 

610 CONTINUE 

PRINT 6, NR, NROUT 

6 FORMAT C2X, 'EOF AFTER N. RECORDS = , ,I10,2X,' NROUT = ' ,110) 
600 CONTINUE 
C 

CALL NTRANS C I0UT2, 1 , 2000 , ROUT , L , 22) 

NROUT = NR0UT+1 
PRINT 5 , NROUT , NR IN 
WRITE C 8 , 400 ) 

WRITEC8 , 5) NROUT, NRIN 

5 FORMAT C2X, 'NUMBER OF RECORDS WRITTEN AND READ' , 2X , 2110) 
CALL NTRANS ( I0UT2 , 3 , 2000 , AOUT , L , 22) 

GO TO 1000 

900 PRINT 350, K , LL , NR, NROUT 
1000 CONTINUE 


FINE DEL MAIN 
STOP 
END 


C 

C ****************-*************^******^^^^********-**^*********** 
C QUEST A SUBROUTINE RIEMPIE IL BUFFER 
SUBROUTINE FILBUFCLL , I 0UT2) 

C 

COMMON /DAT I /AC SOI ), AOUT ( 100) , NROUT 
CHARACTER*80 AOUT 


C 


DOUBLE PRECISION A 


LL=LL+ 1 



n n n n 


IFCLL.LT. 100) RETURN 
COLL NTRONS ( I0UT2, 1 , 2000 , OOUT , L , 22) 
IF(L.LT.O) PRINT 10, NROUT 
NROUT = NROUT + '1 
NO=CNROUT/50)*SO-NROUT 
IFCNO.EQ.O) PRINT 20 , NROUT 
10 FORMAT C 2X , ' HARDWARE ERROR OT RECORD =',110) 

20 FORMAT ( 2X, 1 NUMBER OF RECORDS WRITTEN =',110) 

LL= 1 


RITORNO ML MO IN 
RETURN 

FINE DI FILBUF 
END 



FTN7X , 5 

PROGRAM DEC00EREGRE5SC ) , < SQ0323 . '! 415 > 

C 

C xx^xxx**-**-*x**************-**-k^**-*******-x*-x****.*-**-*.*.***-x.*****-x******* 


C * DECOD1FICP SEQUENZIALE DEI DATI REGRE55 DO NA5TR0 * 
C XXXXXXXXXXXXXXXXX*XXXXXXXXXX*-X-XXXXXXXXXXXXXXX*XXXXXXXXXXXXXX-XX-XXXXXX 
C * * 
C * QUE5T0 PROGROMMA OPERA 5UI DATI REGRESS. * 
C * * 
L * ESSO LEGGE DA NOSTRO RELORDS DI 4000 PAROLE, LI DECQDIFICA * 
C * ED IMMAGAZZINO 5U NOSTRO D'USCITO COME INTEGER *4 . * 
C * * 


C x-x*x-x**xxxx*xx-x*x*x**xx-x*xxxx*-x-x*-x-x'x-x**-x***-xxxxxxx-xxxx*xxx********** 

c 

C DICHIOROZ I ONI 

DIMENSION I5KC40, 100) 

COMMON I SKI C240) , IP, IRM, IRC, IC3, IPRINT1 , IPRINT2 , I TAPE I 
COMMON IWR1 ,NIN, ILM1 , ILM2 
INTEGER*4 LMATC160.15) 

DIMENSION IAC4) 

C 

C E5EGUIBILE 

C I NFORMOZ I ON I DI INGRE55Q/U5CIT0 
111 WRITE Cl, 11) 

11 FORMAT C IX, 21 / ) 

WRITEC1,*) " 

WRITEC1,*) ' DECODIFICO SEQUENZIALE DEI DATI REGRESS’ 

WR I TE C 1 , X ) ' XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX-XXX-X-X ' 

WRITEC1 ,x) 

WRITEC1 ,*) 

WRITEC 1 , x) ' QUE5T0 PROGROMMA LEGGE DO NOSTRO RECORDS DI ' 

WRITE C 1 , x) ' 4000 PAROLE FORNENDO IN USCITO, 0 5CELTQ : ' 

WRITEC 1 , x ) 

WRITE ( 1 , x) ' 1) 5EMPLICE LISTOTO DEI DATI' 

WRITEC 1 , x) ' 2) LISTOTO DEI DATI DECODIFICATI ' 

WRITE C 1 , x ) ’ 3) IMMPGPZZI NQMENTO DEI DATI DECODIFICATI' 

WRITEC1 ,12) 

12 FORMAT C IX , / ,80'x' , / ) 

WRITEC 1 ,' C "RECORD REGRESS INIZIOLE: " ) ' ) 

READC 1 ,x) NIN 

WRITEC 1 ,’ C “RECORD REGRESS FINALE: ")') 

READC 1 , x) NFI 

WRITE Cl , ' C "UN I TO' ' NOSTRO IN LETTURO: ")') 

READ Cl,*) I TOPE 
REWIND ITOPE 

WRITEC 1 , ' C "VUOI IL SEMPLICE LISTOTO? <1=5I>: ")') 

READC 1 , x) IPRINT 
IF C IPRINT . EQ . 1 ) THEN 
WRITE Cl , ' C" UNI TO" 5T0MPANTE? : ")') 

REODC 1 ,x) IWR 
WRITE C I WR , x ) 

WRITEC1 , ' C "DE5IDERI SOLO QUESTO STQMPO? <1=SI>: " ) ’ ) 

REODC 1,*) IONLY 



n n 


ENDIF 

I WRITE(1 ',^C "VUOI^MMOGOZZINQRE I DOTI 5T0ND0RD? <1=SI): ")') 

REPDC1,*) IPRINT2 
TFCIPRINT2 EQ.1) THEN 

WRITE Cl , ' (“UNITO" NOSTRO IN SCRITTURO?: 11 ) ' ) 

REODC'I,*) I TOPE I 

^RITECl]*?’ 'OTTENToI^HPI FGRNITO LO STLSSO UNITO" NOSTRO' 
WRITEC1,*) 'SIR IN 5LRITTUR0 LHE IN LETTURO ! 

GO TO 111 
ENDIF 

WRITEC1 , ' ("VUOI UNO STOMPO DEI DOTI STONDORD? <1=SI>: ">') 

REODC 1 , *) I PR INTI 
IFCIPRINT1 .EQ.1) THEN 
WRITEC 1 , ' C "UNITO' ' bTRMPRNTE:? : ")') 

REODC'I, *) IWR1 
WRITEC IWR'l ,*) 

ENDIF 

ENDIF 

RICERCO DEL RECORD INIZIOLE RICHIE5T0 

IFCNRC N NE.O) WRITEC1,*) 'RICERCO DEL RECORD INIZIOLE' 

DO 1 1=1 , NRC 

COLL XTOPEC1 , I TOPE , I5K , 4000 , LEN ) 

I F ( LEN . NE . 4000 ) THEN 

100 F0RM0TC2X?°ERR0RE in LETTURO REC# ',15,' LUN . # M5) 

ENDIF 

1 CONTINUE 

200 FORMAT C ' NUMERO DI RECORDS REGRESS 50LT0TI : ' , 15) 

WRITEC 1,300) 

300 FORMAT (IX, / , / ,80'*' ) 

PER CIOSCUN RECORD RICHIE5T0 
DO 2 I I =NIN , NFI 
IRC= IRC+1 

LETTURO DI UN RECORD 

COLL XT0PEC1 , I TOPE, ISK , 4000 , LEN) 

IF (LEN . NE . 4000) THEN 
IF (LEN.EQ.-1) THEN 
WRITEC 1 ,*) 'EOF IN LETTURO' 

GO TO 10 

WRITE Cl, 100) 'ERRORE IN LETTURO RECORD #',II, 

GO TO 2 

WRITEC 1,*) 'LETTO DOLLfl UNITO' ", ITOPE , ' IL RECORD #',II 


LUN . = ' , LEN 
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C 5TRMPR IL RECORD REGRESS SE RICHIE5T0 
IFCIPRINT.EQ.1) THEN 

IFCIWR.NE.1) WRITEC 1 ,*■) '5CRITTURR SULLR UNITR ' * ' , IWR , 

*' REC tf , II 

WRITEC IWR, 400) II 

400 FORMRT ( 1 X , ' RECORD # ',15) 

WRITEC IWR, 500) C ( ISKC I , J) , I =1 , 40) , J = 1 , 100) 

500 FORMRT C40R2) 

WRITEC IWR, 600) 

600 FORMRT C / ) 

IF C I ONLY . EQ . 1 ) GO TO 2 
END IF 

PER CIR5CUNR COLONNR DEL RECORD 
DO 3 I CL = 1 ,100 

CONTROLLO SUGLI I DENT I F I CRTOR I 

CRLL SMOUE Cl SKC 1,1 CL), 1,8,10,1) 

DECODECS, 700, IR,ERR-20) IC1 , IC3 
700 FORMRTC 15, 13) 

IFCIC1 .NE.IC2) THEN 
IC2=IC1 
IP=1 
ENDIF 

TRRNSFERIMENTO DI UNR COLONNR DI 'ISK' IN ’ ISK1 ' CONTROLLRNDO , 

DRTO PER DRTQ , SE IL BUFFER E' PIENO CTEST SU IP) 

CRLL SMOVE CI5KC1 , I CL), 3, 32, I SKI ( 1 ) , I P ) 

CRLL CONTROLLO CLMRT) 

CRLL 5M0UE C ISK C 1 , I CL ) , 33 , 56 , I SKI C 1 ) , IP ) 

CRLL CONTROLLO CLMRT) 

CRLL 5M0VEC ISK Cl , I CL) ,57,80, I SKI ( 1 ) , IP) 

CRLL CONTROLLO CLMRT) 

SI RIPETE PER UN' RLTRR COLONNR 
3 CONTINUE 

51 RIPETE PER UN RLTRO RECORD 
2 CONTINUE 

OPERRZIONI FI NOLI 
10 REWIND ITRPE 

IFC I ONLY . NE . 1 ) THEN 

SI I MMRGRZZ I NRNO I DOT I RIMRSTI IN MEMORIR 
IFC IPRINT2.EQ. 1 ) THEN 
I LM2= ILM2 4- 1 

CRLL XTRPE C 2 , I TRPE 1 , LMflT , 4800 , LEN ) 

IFC LEN . NE . 4800 ) THEN 

WRITEC1 ,*) 'ERRORE IN SCRITTURfi SU ',ITRPE1 
WRITEC 1 , ■*) ' LMRT #',ILM2,' LEN =',LEN 
ENDIF 

WRITE Cl,*) '**SCRITTO SU' , I TRPE 1 , ’ #', ILM2 ,' RECORDS STRNDRRD ' 



ENDIF 

IF ( IPRINT1 . EQ . 1 ) THEN 
I LM1 = ILM1 + 1 

WRITE Cl ,*) 'SCRITTURO SU' , IWR1 , ’ DEL RELORD #',ILM1 
DO 4 MM= 1 , '1 60 

4 WRITEC IWR'l ,800) MM , ( LMOT ( MM , NN ) , NN = 1 , 15 ) ^ 

800 FORMOT CI3, 1 6, 17, 19, 17, 19, 18, 19, 1 6, 17, 18, 517) 

ENDIF 

C 

f SCRITTURO DELL' EOF 5UL NOSTRO IN U5CIT0 
IFC IPRINT2.EQ. 1 ) THEN 

WRITEC 1 ,' C "UUOI UN OLTRO NOSTRO RtbRESb? C1=bl): ")') 

RERDC1,.*) NY 

IF C NY .EQ. 1 ) GO TO 111 

ENDF ILE CUNIT= I TRPE1 , IOSTOT=K<K , ERR=30) 

WRITEC'I , *) '**5LRITT0 EOF SL'.ITOPEl 
REWIND IT0PE1 
ENDIF 
ENDIF 


C FINE REGOLORE DELL' 0N0LI5I 
WRITEC'I ,300) 

WRITEC'I,*) '*****' 

WRITEC1,*) ' ***LETT I DOMTOPE,’ ’ , NFI -NIN + 1 , ' R. RECORDS' 
IFC IPRINT .EQ. 1 ) THEN 
WRITEC 1,*) '*****' 

WRITEC'I,*) ' 5CRITTI SU'.IWR,' ',IRL,' R. RECORDS' 

ENDIF 

IFC IPRINT2 . EQ . 1 ) THEN 
WRITEC 1 ,*) '*****' 

WRITEC'I,*) ' IMMOGOZZINOTI SU ' , I TflPEI , ' ',ILM2,' 5. RECORDS' 

ENDIF 

IF C IPRINT1 . EQ. 1 ) THEN 
WRITEC 1 ,*) '*****' 

WRITEC 1,*) 'SCRITTI SU' , IWR1 , 

ENDIF 

WRITEC'I,*) '*****' 

WRITEC 6 , 300) 

WRITEC 6,*) '*****' 

WRITE C6 , *) ' ***LETTI DO',ITOPE 
IF C IPRINT . EQ. 1 ) THEN 
WRITEC 6,*) '*****' 

WRITEC6,*) 'SCRITTI SU',IWR,' 

ENDIF 

IFC IPRINT2 . EQ . 1 ) THEN 
WRITEC6,*) '*****' 

WRITEC 6,*) ' IMMObOzZINOTI SU 

ENDIF 

IFC IPRINT1 .EQ. 1 ) THEN 
WRITEC 6,*) '*****' 

WRITE C 6 , *) 'SCRITTI SU'.IWRI, 

ENDIF 

WRITEC 6,*) '*****' 


' , ILM'1 , ‘ 5. RECORDS' 

’ ' ,NFI-NIN + 'l , ' R. RECORDS' 

' , IRC, ' R. RECORDS' 

I TOPE 1,' ' , I LM2 , ' 5. RECORDS' 

' , ILM1 , ' S. RECORDS' 
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WR I Tt ( I , yoo ) 

900 FORMOTC IX, / , 50 ' * * , ’FINE REGOLPRE DELL' ' ONRL. ISI ’ ) 
REWIND I TOPE 
REWIND I TOPE 1 
STOP 
C 

C ERRORS D I DECQDIFICR 

20 WRITEC I , *) ' **ERRORRE DI DECODIFICR** ’ 

STOP 


ERRORS IN EOF 

30 WRITE Cl,*) ' ERRORS SLR, I TREE!,’ IN EOF tt ' , KKK 
STOP 

FINE DEL MR IN 
END 


QUESTR SUBROUTINE SLITTR L ' INDICE DI P05IZI0NE NEL BUFFER 
E LO DECODIFIER 5E QUESTO E' PIENO 
SUBROUTINE CONTROLLO ( LMRT ) 


DICHIRRRZIONI 

COMMON I SKI (240) , IP, IRM, IRC, IC3, IPRINT'I , IPRINT2 , I TAPE"! 

COMMON IWR1 ,NIN, ILM1 , ILM2 
INTEGER*^ LMRT (160,15) 

INTEGERS LYD , LHM5 , LTRM , LE , LB , LFRQ 

DOUBLE PREC IS ION DT I M , DOBS , DFRQ , DPR5 , DCOM , DRTL , DHRG , DDEC , DRZM 
DOUBLE PRECISION DELU , DG 1 , DGM3 , DGM'l , DRE5 , DRJC , DWGH , DCRS , DMOD 
CHRRRCTER 01,02,03,04,0(10) 

SLITTRMENTO DELL' INDICE DI POSIZIONE NEL BUFFER 
IP=IP+24 

CONTROLLO SUL RIEMPIMENTO DEL BUFFER 
IFCIP.EQ.481) THEN 

IL. BUFFER E' PIENO: CONTIENE TUTTI I DAT I RELRTIUI RD UNO MI 51 IRR 
51 EFFETTUR LR DECOD I FI CR 
IP=1 

DECODE (480, 100, ISK'l ,ERR=900) DTIM,01 , LYD , LHMS , Cfl(K) ,K= ! , 10) , 
+IB, IT, IR, I M, 02, 03, 04, DOBS, DFRQ, DPR5, DCOM, DRTL, DHRG, DDEC, DRZM, 

+ DELU , DG1 , DGM3, DGM'l , DRES , DRJC , DWGH , DCRS , DMOD 
100 FORMAT (D24 .18,02,15, 17, 1002, 412, 302 , 17D24. 18) 

CONTROLLO SULLR BON TO ' DEI DRTI 

I F C DRJC . NE . 0 . DO . OR . DWGH . NE . 0 . DO ) THEN 
ICO= IRC fNIN- 1 

WRITEC 1 , *) ' RECORD ICO, IC3 

WRITE C 1 , *) ' ERRORE NEI DRTI: REJCOD =',DRJC,' WEIGHT = ' , DWGH 

WRITEC1,*) ' ' 

END IF 



C ELOBOROZIONE DEI DOT I 
LHM5=LHM5/ 10 - 
LTRM=IT*'100Q0. + IR* I00 1 +IM 
COLL ELEVCIT.DHRG.DDEL.EL) 
IELUF=D61* 100 . 

I ELDW=EL*1 00 . 

LE=IELUP* 10000. +IELDW 
LB=DRE5*10OOOO0 . 

I_KRQ=DFRQ 


f TROSFFRIMENTO DEI DOT I NEL BUFFER DI UbClTO 

IFCLMOTCJRM 1 ) . EO . l.YD . OND . LiiPT C IRM , 2 ) . EU . L.HM5 
+ .PQ. i .FRQ . OND . LMOT C IRM , 4 ) . EQ . LTRM . OND . LMOT C IRM 

C MISURO CONTEMPORPNEO OLLfl PRECEDENTS 
IF( IB . EQ .11) THEN 

C 5 -BOND 

LMOT C I RM , 6) =LB 
l_MOT C IRM , 9 ) = DOBS 
ELStr. 


MND . LMP i ( I RM , 3 ) 
5i.EQ.LE) THEN 


C X-BOND 

LMOT ( IRM , 7) =LB 
LMOT (IRM, 10)=D0B5 
IFCIB.NE.21) THEN 
L MOT C IRM, 12) =939333999 
ELSE 

TBD=LMOT ( IRM , 6) -3 . /II . *LMOT( IRM , 7) 

LMOTC IRM , 8 ) =TBD 
END IF 
END IF 

C TROSFERIMENTO IN UNO NUOVP RIGO DEL BUFFER LMOT 
IRM= IRM+1 


C 

c 


CONTROLLO SUL RIEMPIMENTO DEL BUFFER DI USCIiO 
IFCIRM.EQ.1B1) THEN 
IF( IPRINT2.EQ. 1 ) THEN 


ILM2 = I LM2 + 1 

COLL XTOPE C 2 , I TOPE 1 , LMOT , 4800 , LEN ) 

IFCLEN . NE . 4800) THEN 

WRITEC1 ,*) 'ERRORE IN bCRITTURO bU'.ITOPEI 
WRITE Cl,*) 'LMOT #',ILM2,' LUN = ' , LEN 

WRITEC1,*) ' **SCRI TTO SU' , I TOPE 1 , ' 5. RECORD 


tr 


ILM2 


I 


>00 


ENDIE 

I F ( I PR I NT1 . EQ . 1 ) THEN 
ILM1 = ILM1 +• i 

WRITEC 1 ,*) 'SCRITTURR bU',IWR1, 


5. RECORD tf , ILM1 


DO 1 MM = 1,1 bO % 

WRITEC IWR1 , 200 ) MM , ( L MOT ( MM , NN ) , NN- 1,15) 
FORMOTC 13, 16, 17, 13 ,17, IS, 18,13,16, 17, 18, 517) 


ENDIE 

DO 2 MM = 1 , 160 
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DO 2 NN = 1,15 
2 LMOT C MM , NN ) =0 

IRM= 1 
ENDIF 

C 

C MI5UR0 NON CONTEMPORONEO OLLO PRECEDENTS 
LMOTCIRM, I ) = LYD 
LMOT ( IRM ,2) =I.HM5 
LMOT C IRM, 3) = DFRQ 
LMOTC IRM, 4)=I_TRM 
LMOT ( IRM , 5)=!_E 
LMOT C I RM , 1 I ) = DR JC*1 0 . **6 
LMOT C I RM , 1 2 ) = OWGH* '1 0 . **6 
I F C I B . EQ . 1 1 ) THEN 

E S-BOND 

LMOT C IRM , 6 ) =LB 
LMOT C IRM, S)=DQBS 
ELSE 


C X-BOND 


LMOT C IRM , 7 ) =LB 
LMOT C IRM ,10) =D0B5 
ENDIF 
ENDIF 


RITORNO OL MO IN 
ENDIF 
RETURN 

ERRORE DI DECODIFICO DEL EUFFER 

200 WRITE ( '1 , *) '**ERRQRE DI DECODIFICO DEL RECORD**' 
STOP 


FINE DI CONTROLLCJ 
END 


***********************************************************#**#****- 5 * 


QUESTO ROUTINE COLCOLO L' EL.EUOZIONE 
SUBROUT I NE ELEV (IT, DHRG , DDEC , EL ) 


D I CH I OROZ I ON I 

DOUBLE PRECISION DHRG, DDEC 


RICERCO DELL.O LOTITUDINE DELL.0 STOZIONE IN OUFSTIONF 
IF ( I T . EQ. 12) THEN 
L0TD=35 
LOTP= 17 

ELSE 

I F C I T . EQ .14) THEN 
L0TD=35 
L0TP=2S 

ELSE 

IFC IT .EQ. 42) THEN 
LOTD= -35 
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LPTP = 2^4 


ELSE 

IFfIT.EQ.43) THEN 
LPTD = -35 
LRTP=24 

ELSE 

IFCIT.EQ.61) THEN 
LRTD=40 
LllTP = 25 

ELSE 

IFCIT.EQ.63) THEN 
LRTD=40 
I.PTP = 25 
END IF 
END IF 
ENDIF 
ENDIF 
ENDIF 
ENDIF 
C 

C CRLCGLR L ' ELEVOZIONF 
PIG=3 . 141592654 
CON=PI6/iaO - 
IFCLRTD.GT.O) THEN 
PLAT =CGN*(LRTD+LRTP / SO . ) 

ELSE 

RLRT=C0N*CLRTD-LRTP/60 . ) 

ENDIF 

RHRG=DHRG*CON 
RDEC = DDEOCGN 

5 1 NEL = COS ( RDEC ) *C05 C PLOT ) *LU5 RHRG ) +5 1 N C RDEC ) *5 1 N ( PLOT ) 
EL=RSIN(SINEL) 

EL=EL/CON 

RITORNQ PL MR IN 
RETURN 

FINE. DI ELEV 
END 
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<860323. 1415) 



DICHIRRAZIONI 

I NTEGER*4 LMRTC160, 15) , I YD, ZERO 
DIMEN5I0N NRTR (48,24) 

CHARACTER B 

INFORMAZIONI DI I N6RE5S0 / U5C I TP 
WRITE Cl ,100) 

100 FORMAT C 1 X , 20 / , 80 ' * ' ) 

WRITEC1,*) ' LISTR RIPS5UNTIVA DEI DRTI STRNDRRD ' 

WRITE! 1 !*) 

WRITE!'!,*) ' QUESTO PROGRRMMR PRODUCE ITERRTIURMENTE , ' 
WRITE!'!,*) ' PER 06NI GIORNO DI DRTI STRNDRRD RICHIE5T0,' 
WRITE! 1,*) ' UNR LISTR RIASSUNTIUR DEL NUMERO DI DRTI' 
WRITE!'!,*) ' SEPRRRNDOLI PER ORR, 5TRZI0NE , BANDA E MODO . ' 
WRITE! 1 ,101) 

101 FORMRT! IX, / ,80'*' , / ) 

10 WRITE! 1 ,'! "NRSTRO DI LETTURR CENTER 8 OR 9) #: " ) ' ) 

RERD Cl,*) I TAPE 

IF! I TAPE . NE . 8 . RND . ITRPE . NE . 3) GO TO 10 
REWIND ITRPE 

30 WRITE! 1 ,' C "GIORNO INIZIRLE CDR 1 TO 365): " ) ' ) 

RERD! 1 , *) I DRY I 

IFCIDAYI .LE.O.OR. IDRYI .GE.366) GO TO 30 
40 WRITE! 1 ,' C "GIORNO FINALE CDR 1 TO 365) #: H ) ' ) 

READ!'],*) IDRYF 

IF! IDRYF. LE.O.OR. IDAYF. GE.366) GO TO 40 
50 WRITE Cl , ' C "ANNO: 19 ")') 

RERD Cl,*) I YEAR 

IF! I YEAR .GE . 100) GO TO 50 

WRITE Cl , ' C "STRMPANTE ft: ")') 

RERD! 1 , *) IWR 
ZER0=0 

PER CIRSCUN GIORNO RICHIE5T0 
DO 8 I D= IDRYI , IDRYF 
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IYD=IYEAR*1000. *ID 

I N I Z I RL I ZZOZ I ONE 
NRC = 0 

DO 1 11=1, 24 
DO 1 JJ=1,4fl 

1 NRTRCJJ, II )=0 

CALL EXECC3, ITRPE+0200B) 

I PAGE = 12 

WRITFC IWR , 200) IPRbE 
200 EORhRT C 02 ) 

LFTTURR DI UN RECORD 

2 COLL. XTRPEC 1 , I TOPE , LMOT , 4800 , L.EN ) 

NRC = NRO 1 

IF(LEN.EQ.-"n THEN 
WRITEC1,*) 'EOF IN LETTURO’ 

GO TO 60 
END IF 

IFCLEN.NE.4800) WRITE Cl , *0 ' **ERRORE IN LETTURO: LUN . = ' , LEN 

PER EI05EUN0 RIGA DEL RECORD 
DO 3 1=1,160 
I RG= I 


RICERCO DEL GIORNO RILHIESTO 
I GO=LMRT (1,1) - I YD 
IFCIGO.LT. 0) GO TO 3 
IFC IGO.GT.O) GO TO 60 

0NQLI5I DEI DOT I PER LO MTRICE 

I RM=LMQT C I , 4 ) -LMOT C I , 4 )/ 1 0000*1 0000 
I R = I RM / 1 00 
I M= I RM - I R*1 UO -10 
IFCIR.EQ.12) THEN 
I RR = 1 

ELSE 

IFCIR.EQ.14) THEN 
IRR = 2 

ELSE 

IFCIR.E0.42) THEN 
I RR = 3 

ELSE 

IFCIR.EQ.43) THEN 
IRR = 4 

ELSE 

IFCIR.EQ.61) THEN 
I RR = 5 

ELSE 

IFC IR . EQ. 63) IRR=6 
ENDIF 
END IF 
ENDIF 
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ENDIF 
END IF 

ih=lmpt c i , 2 ) / lonon 

I MS = LMPT C I ,2)-IH*10000 
I B T = I H*SO +• I MS / 1 00 + 1 

c 

c riemhimento dellp listp 

12= C IBT-1 ) / 60 . + I 
IB=1 

I UCIRR-1 IM-1)*2+IB 

JFCLMOTC I ,6) . NE. ZERO) NRTRC I 1 , I2)=NRTRC I I , 12) + I 
I B = 2 

I 1=CIRR- 1 )*8f C IM- 1 )*2+IF) 

IFCLMOTCI ,7) . NE . ZERO) NRTRC 1 1 , I2)=NRTRC I I , 1 2 ) + I 

51 R I PETE PER UN'RLTRR RIGA 

3 CONTINUE 

51 RIFETE PER UN RLTRG RECORD 
GO TO 2 

FINE DELLO RICERCR DEI DRTI 
60 IYERR1=LMPTC IRG, 1 )/1000. 

I DRY 1 =LMRT ( IRG , 1 ) - I YERR1*1 000 . 

WRITEC 1 ,*) 'FINE RICERCR IL',IDRY1,' DELL*' RNNO ' , I YERR'I 

5CRI TTURR DELLO MRTRICE DI U5CITR 
DO 4 JJ=1 , 48 
DO 4 11= 1 , 24 

IFCNRTRCJJ, II) .NE.O) IPRINT1=1 

4 CONTINUE 
IFC IPRINT'I . EQ. 1 ) THEN 

I PRINT 1 =0 
WRITEC IWR, 300) 

300 F0RMRTC2X, 132' ) 

WRITEC IWR, 400) IYEOR,ID 

400 F0RMPTC2X, ' YEOR = ',I3,5X,'DPY = ' , I 4 , 3X , ' HOUR = 012 

5 : 6783 10 11 : 12 13 14 15 16 17 : 18' 

13 20 21 22 23 NS DRY',/,/) 

DO 7 JJ=1 , 48 
DO 5 11=1,24 

5 IFCNRTRCJJ, II ) .NE.O) IPRINT=1 
IFC I PRINT . EQ. 1 ) THEN 

IPRINT =0 

is = c jj - 1 ) /a. t-i 

IFC IS . EQ . 1 ) THEN 
155=12 

ELSE 

IFC 15. EQ. 2) THEN 
155=14 

ELSE 

IFC IS . EQ . 3) THEN 
I SS = 42 
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ELSE 

IF C 15 . EQ 

. 4 ) THEN 

155 = 43 

ELSE 

IFCIS.EQ.5) THEN 

155=61 

ELSE 

IF C 15 . 
ENDIF 

EQ.S) 155 

FNDIF 

ENDIF 


ENDIF 
ENDIF 
IP= JJ/2 . 

IFC IP*2.EQ. 

,U) THEN 

B= ’ X ’ 
IBB =2 

ELSE 

B= ' S’ 
I BB = 1 



END IF 

IMM=C JJ-IBB-CIS-1 )*8)/2+1 


I TOT = 0 


DO 6 LL=1 , 24 

6 I TOT = ITOT +NQTRC JJ , LL ) 

IFCIS'l .NE. 15) WRITECIWR,*) 

151=15 

WRITEC IWR , 500) 1 55,8, I MM, CNHTR(JJ,LL) , LL-'i , 
500 F0RMPTC2X, ' DOP : D55 =',I3,' Band-Mode 

: ' , 613 , ’ : ', 613 ,' : ', 613 , 18 ) 


24) , ITOT 

' , 01 , 12 


ENDIF 

7 CONTINUE 

WRITE ( IWR ,#) 

WR I TE C I WR , 300 ) 
ENDIF 


SI RIPETE PER UN RLTRO 6IQRNQ 
IFCLEN.E0.-1) GO TO 70 
8 CONTINUE 


FINE RE60L0RE DELL' 0N0LI5I 
70 WR I TE C I WR , 200 ) IPP6E 
WRITEC IWR , 200) IPOGE 

WRITEC1,*) . 

WRITEC'l,*) '**FINE REbGLOF^E DELL' PNPLIbI** 

STOP 


FINE DEL MO IN 
END 


, 4X , 613 , 
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‘N7X , 5 


PROGRAM PL07REGRE55 C ) , 


<860323. 14'I5) 


* TROCC I OMEN TO SEQUENZIALE DEI DOT I STANDARD . ^ ^ . * 

* * 

* 01 IE 5 TO FROGRQMMP OPERA EOT HAT I STANDARD * 


* ESSO TRACE JO TIJTTI I DAT I RELOT K 1 1 01... GIORNO, ML LA STOZIONE ED * 

* PL MODO RICHIE5TU . * 

* QUEST I VENGONO TROCl I OT I SEPAROTI PER BANDA IN. THE GRGFICI . * 

* NEL MOSS J MO FORMOTO. * 






DICHJORAZIONI 

EMO B5 C 720 ) , BX ( 720 ) , BD C 720 ) , BT C 720 ) 

INTEGER-x-4 LMPTC 160, IS) , I YD 
DIMENSION VC 23) , 1ST (S ) , I MODE (4 ) 

CHARQCTER*50 N1,N2 
INTEGER YN6 
CHPRPCTER*2 ANNO 

E5EGUIBILE 

INFORMOZIONI PER INGRESSO/USCITO 
WRITE Cl, 100) 

1 00 FORMAT C 1 X , 20 / , 60 ' * ' , / ) 

WRITEC'I,*) ' TROCCIAMENTO SEQUENZIALE DEI DATI STANDARD' 

WRITE C 1 !*) 

WRITE Cl , *) 

WRITEC'I,*) ' QUESTO PROGROMMA TROCC 10 TUTTI I DATI' 
WRITEC'I,*) ' RELOTIVI A CIOSCUN GIORNO RICHIESTO' 
WRITEC'I,*) ' IN SUCCESS I ONE PER STOZIONI E MODI' 

WRITEC'I ,*) 

WRITEC1,*) ' I GROFICI SONO 3 PER OGNI FOGLIO' 

WRITEC I,*) ' II. FORMOTO E'* IL MASSIMO' 

WRITEC'I ,*) 

WRITEC'I,*) ' ' 

WRITEC 1,*) ' PENNINI DEL. PLOTTER: ' 

WRITEC'I,*) ' 1: PER LE I NTE5T0ZI0NI ' 

WRITEC'I,*) ' 2: COMPLESSO 10' 

WRITEC'I,*) ’ 3: COMPLESSO 40' 

WRITEC'I,*) ' 4: COMPLESSO 60' 

WRITEC 1 ,*) 

WRITE C 1 , 101 ) 

101 FORMAT C 1 X , / ,60' *' , / ) 

10 WRITEC 1 ,' C "NOSTRO DI LETTURO C8 0 S) #: ")') 

READ C'l,*) I TOPE 

IFC I TOPE . NE . 8 . AND . I TOPE . NE . 3 ) GO TO 10 
REWIND I TAPE 



30 WRITEC 1 , ' C "GIORNQ INIZ.IPLE (DPI 
REP DC 1 , *) I DRV I 

I F ( I DRY I . LE . 0 . OR . I DRY I . GE . 36b ; 

40 WRITEC !,' ( " GIORNO FINALE (DR ' 

RERDC I , *) IDRYF 
IF C I DRYF . I..E . 0 . OR . IDRYF . bE .36b ■ 

50 WR I TE C 1 , ’ C " RNNO : 1 3 " ) ' ) 

RERDC I,*) I YEAR 
IFC I YEAR . GE . 100) bO TO 50 
WRITEC RNNO, 102) IYERR 
10? FORMAT (12) _ , 

WRITEC 1 , ' C “ VUOI L” ECO SULLA b? <Y=!>: 

RERDC!,*) YN6 
IFCYN6.EQ. I ) THEN 
IPRGE = "12 

WRITEC6,200) IPRbE 
200 F0RMRTCR2) 

ENDIF 
NR = 0 


1 R 365) X: " ) ' ) 

GO TO 30 
R 365) X: " ) ' ) 

GO TO 40 


NM = 0 

DO 1=1,6 

NR = NR + 1 . 

WRITEC1,*) ‘ ***** 0 PER TERM^NPRt: 
WRITE (1 , ' ( "STPZIONE RICEVENTE? <12 
RERDC 1 , *) ISTCNR) 

IFCI5TCNR) .EQ.O) THEN 
NR=NR-1 
GO TO 60 

ENDIF . 

I F C 1 5T C NR ) . NE . 1 2 . RND . 1ST CNR) . NE . 1 ^ 

+ ISTCNR) . NE. 43. RND. ISTCNR) .NE. 61 .RND. IS 
NR=NR-1 


14,42,43,61,63): ")') 


RND . I ST C NR ) . NE . 42 . RND 
CNR) . NE . 63) THEN 


WRITEC 1 
WRITEC1 
ENDIF 
END DO 
60 DO I = 1 , 4 
nm=nm+i 

WRITEC! , 
WRITEC 1 , 


,*) 

,*) 


- *****5TRZI0NE NON RICONUSCIUTR !*** 


*) ****** 0 PER TERM IN ARE ’ 

' C "MODO RICHIE5T0? < DR 1 fl 
RERDC 1,*) IMODECNM) 

IFC I MODE CNM) . EQ.O) THEN 


4 > : " ) ' ) 


NM=NM-1 
GO TO 70 
ENDIF 

IFC IMODECNM) 
NM=NM- 1 
WRITEC 1 ,*) 
WRITE Cl ,*) 
ENDIF 
END DO 
70 CONTINUE 


LT . 1 .OR . IMODECNM) . GT .4) THEN 
*****MQDO NON R I C0N05C I UTO ! ***** ' 


C 
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r PER CIR5CUN GIORNO RICHIFbTO 
DO 7 T D= I DRY I , IDRYF 
I YD= I YEPR* TOGO . +ID 


PER CIRSCUN MODO RICHIEbTO 
DO 5 IMM=1,NM 
IM= IMDnFC I MM) 

PER CJRSCUNR STRZIGNF RICHIES7R 
DO 6 K=!,NR 
1 5= I ST ( K' ) 

INIZIRL I ZZRZIONF 
IP0 = 0 
NRC = 0 

DO 1 I = ’1,720 
BSC I ) =0 . 

BX C I ) =0 . 

BDC I ) =0 . 

BTC I ) = 0. 

1 CONTINUE 
IFCYNS.EQ. 1 ) THEN 

WRITEC6,*) 

WRITE C 6 , -*) 

WRITEC6,*) ' DOY =',ID,' RNNIJ =',IYERR, J bTRZIONE RIC.=',I5, 
+ ' MODO = ’ , I M 
WRITE CS,*) 

WRITEC6 , #) ' MINUTO S-BRND X-BRND NON RIPORTRTE ' 

END IF 

WRITEC 1 ,*) 'DOY = ' , ID , ' RNNO =',IYERR,' STRZIONE RIC.=',IS, 

+ ' MODO = ’ , IM 
REWIND ITRPE 

LEGGE UN RECORD 

2 CRLL XTRPEC1 , ITRPE , LMRT , 48U0,LEN) 

NRC = NRC +1 

IFCLEN.NE.4800) THEN 
IFCLEN.EQ.-1) THEN 
WRITE C l , * ) ' -**EOF SU NOSTRO *, ITRPE 
GO TO 4 

ELSE 

WRITE Cl,*) '**ERRORE IN LETTURO: LUN= ' , LEN , ' SU ITRPE 
END IF 
END IF 

PER TUTTE LE RIGHE DEL BUFFER 
DO 3 1=1,160 
IRG= I 


RICERCR DEL GIORNO RICHIE5T0 
I GO = LMRT C I , D-IYD 
IFCIGO.LT.O) GO TO 3 
IFCIGO.GT.O) GO TO 4 



c CONTROLLO DHL MODO E DELLfl STRZIONE 

IRM = I.MRT( I , 4) -LMRT (1,4)/ 1 OOOC*1 OOOO 
IR= IRM/ 100 
IM'I = IRM- I R* 100- 10 
C 

C SE IL DP TO E* QUELLO LERLRTO 

IFC IM.EQ. I Ml . PND . Ib.EQ. IR) TEEN 

C 

f IMMPGPZZINRMENTO DEL DPTO PER LR GRRFILP 
IP0=IP0+1 

BSC IPO ) =LMPT ( I , hi ) / 1000000. 

IFCBSC IPO) .61. .5) BS ( I PO ) = . E 
I F ( B5 C I PO ) . L T . - . 5 ) BSC IPO) = - .5 
BXC IPO)=LMRT(I ,7)/ 1000000. 

IF(BXCIPO) .LT.-2.S) BXCIP0)=-2.S 
IF(BXCIPO) .6T.2.S) BXC IP0)=£.5 
DD C I PO ) =LMRT C I ,8)/ 1000. 

IFCBDC IPO) .GT .20) BD(IPQ)=2G 
IF C BBC IPO) .LT. -20) BBC IPO) = - 20 
I H=LMRT (1,2) / 10000 
IMS=LMRT(I ,2)-IH*10000 
BT C IPO )=IH*6Q+ IMS/ 100+1 

I F C C BS C I PO ) . EQ . 0 . OR . BX C I PO ) . EQ . 0 ) . PND . YN6 . EQ . 1 ) THEN 
WRITECS, 300) BTC IPO) , LMRT C I ,6) , LMRT C I ,7) 

300 FORMAT C 2X , F5 . 0 , 1 1 1 , 111 ) 

IP0=IP0-1 
END IF 
ENDIF 
E 

C 51 R I PETE PER UN' PLTERR RIGP 

3 CONTINUE 
C 

C 51 RIPETE PER UN RETRO RElORD 
GO TO 2 
C 

C FINE RICERCP DEI DPT I 

4 I YERR1 =LMMT C IRb , 1 ) / lULJO . 

I DRY 1 =LMRT ( IRG , 1 ) - I YEAR 1*1 000 

WRITEC1,*) 'FINE RICERCP IL : ' . I DRY! , ' DELL" RNNO ' , I YERR'l 
C 

C IN P55FNZR DI DPT I NON SI TRRLLIR IL GRRFICO 
IFC BTC 1 ) . EQ . 0) THEN 

WRITE Cl,*) '**NON Cl bONO DP' I PER QUEbTO INbREbbO**' 
IFCYN6.NE.1) GO TO 6 

WRITECS,*) ' **NON Cl 50N0 DR 'I DR GRRF7CRRE PER QUESTO INPUT**' 
GO TO 6 

ELSE 

WRITEC1,*) '**6RRFICI IN IJSCITfl : ' , IPO , ' PUNT I**' 

IFCYN6.E0.1) THEN 

WRITECS,*) ' **GRRFlLl IN UbCI Tfl : ' , IPO , ' PUNT I ' 

ENDIF 

ENDIF 



C TRPCC I PMFNTO GRRFICI 

vc i )=o 

VC2 ) = - .5 

V ( 3 ) = . 5 

vc4>=n 
VCR) =0 

V C 6 .) = 1 440 
VC 7) =50 
VC 8) =30 
VC3 > =350 
V CIO) =50 
VC 1 1 )=2 
VC 12) =3 

v c 1 3 ) =sn . 

V C 1 4 ) = . 02 
VC 1 5 ) = 4 

V C 1 6 ) =5 

V C 1 7 ) = 1 
VC 16) =0 
VC 13) =0 
VC 20) =0 
VC 21 )=□ 

V C 22 ) = 1 
VC23) =3 
VC24 ) =4 

IFC IS.EQ. 12. OR. IS.EQ. 14) VC 25) =2 
IFCI5.EQ.42.0R. I5.EQ.43) VC25)=3 
IFC IS. EG. 61 .OR. IS.EQ. S3) VC25)=4 

V C 26 ) = 1 
VC27 ) = 1 
VC28 ) =0 

V C 23 ) = 1 

WRITEC1 ,*) ' TRRCCIPMENTO PR I MO GRPFICO' 
CPLL SETPRRC V, N1 , N2) 

CRLL EGRRPHY CBT , BS , IP0.1 . ,0, 1) 

WRITE! 1,*) ' TRPCCIPTO PRIMO GRPFICO' 

N2= ' ' 

N 1 = ' X-BRND IHz]' 

VC2)=-2.5 

V C 3 ) =2 . 5 

v c fi ) = i on 

VC 10)= 160 

V C 1 4 ) = . 1 

WRITEC 1 , *) ' TRRCCIPMENTO SFCONDO GRPFICO' 
CPLL SETPPR C V , N1 , N2 ) 

CRLL EGRRPHY C BT , BX , I PO , 1 . ,0,1) 

WRITE Cl,*) 'TRPCCIPTO SECONDO GRPFICO' 

N2= ' ' 

N1 = ' 5-3 / 1 1*X [mHz]' 

VC2 ) = -20 . 

VC3) =20 . 

VC8 ) =1 70 



VC 10) =230 
VC14)=1 

WRITEC1,*) ' TRflCCIPMENTO TERZCi bRWF I LO 
CRLL SETPRRCV.NI ,N2) 

CRLL E5RRPHY ( B l , BD , I PO , I . , 0 , ^ ' _ 

WRI TE C I , *) ' TRRCCIRTO TERZO hRRFIlO' 

C SI RIPETE PER UN' OLTRfl 5TRZI0NE^_ 
IFCYN6.EG. 1 ) WRITE (6, 200) IrRlE 
6 CONTINUE 
C 

C COMPLETRMENTO DFL bRRPlLLJ 

N 2 = ' VOYRGER 19'// WNNO / / ‘ 

N1 = ‘ S-BRND IHr 1 ' 

VC1 )=0 
VC2) = - -5 
VCS) = . 5 
VC4 ) =0 
VC5)=0 
VC6)=1440 
VC 7) =50 

V C 8 ) =30 
VCS) =350 
VC 10) =90 

vei l ) =2 
VC 12) =3 
VC 13) =60 . 

V C 1 4 ) = . 02 
VC 15)=4 
VC 1 6 ) =5 

V C 1 7 ) = 1 
VC 18) =0 

V C 1 9 ) = 1 
VC20) =0 
VC 21 )=1 

V C 22 ) = 1 
VC23) =3 

V C 24 ) =4 

V C 25 ) =0 

V C 26 ) = 1 
VC 27) =1 

V C 28 ) =0 
U C 28 ) = I 

yjRI TE C 1 , *- ) ' COMF1. FTRMENTO DEL GRRFICO’ 
CRLL 5ETPRR C V , N1 , N2 ) 

CRLL E6RPPHY C BT , BS , IP0.1 . ,0,1) 

N2= ' ' 

N1=' X-BRND [Hzl' 

VC2) =-2 . 5 
VC 3) =2 .5 
VCS) =100 

V C 1 0 ) = 1 60 

V C 1 4 ) = . 1 


REbREbn F 
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COLL 5ETP0RCV,N1 ,N2) 

COLL EGROPHY ( BT , BX , IPO, 1 .,0,1) 

N2= ' 1 

N 1 = ' 5 - 3 / 1 1 *X [ mH r 3 ' 

V (2) = -20 . 

0(3) =20 . 

V C 8 ) = 1 70 
GC 10) =230 

V C 1 4 ) = 1 

COLL 5ETPPR CU , N I , N2 ) 

COLL EGROPHY C BT , BD , IPO, 1 . ,0, 1 ) 

COLL ZDLIMCO. ,370. ,0. ,280. , I ERR ) 

IFC IERR.NE.O) WRITE Cl,*) ' ZD! IM ERROR , IFRR 
COLL ZP5PLC370. ,260. ) 

COLL Z V IFWCO . , 370 . , 0 . , 280 . ) 

COLL ZWINDCO. ,370. ,0. ,280. ) 

COLL ZC0LRC1) 

COLL SYMBOL ( 43 . , 25 . , -2 . , ' 00 . 00 ' , 0 . , 5 ) 

COLL SYMBOL 03. , 25 . , -2 . , ' 04 . 00 ’ , 0 . , 5) 


COLL 

SYMBOL Cl 43. 

,25. 

,-2. 

, ' 08 . 00 ’ 

,0. ,5) 

COL.L 

SYMBOL C l 33. 

,25. 

i -2- 

, ' 12.00' 

,0. ,5) 

COLL 

SYMBOL C 243. 

,25. 

, £- . 

,'16. 00 ' 

,0. ,5) 

COLL 

SYMBOL C 293 . 

,25. 

, -2. 

, ' 20 . 00 ' 

,0. ,5) 

COLL 

SYMBOL C 343 . 

,25. 

, -2. 

,' 24.00' 

,0. ,5) 


COLL SYMBOL ( 38 . , 30 . , -2 . , ' - . 5 ' , 0 . , 3 ) 

COLL SYMBOL C 35. , 1 00 . , -2 . , ’ -2 . 5 ' , 0 . , 4 ) 

COL.L SYMBOL C 38. , 170. , -2. , ' -20. ' ,0. ,4) 

FDOY = FLOOT C ID) 

FDSS= FLOAT CIS) 

FMOD = FLOOT C IM) 

COLL SYMB0LC70. ,260. , -3. , ' DOY =',0.,5) 

COLL NUMBER 05 . ,260. , -4. ,FD0Y,0. ,3, ' CI3) ' ) 
IFCNR.EQ.1) THEN 

COLL NUMBER C 155 . , 260 . , -4. ,FD55,0. ,2, ' (12) ' ) 
COLL SYMBOL Cl 30. ,260. , -3. , ' DSS 5) 

END IF 

COLL 5YMB0LC190. ,260. ,-3. , 'MODE =’,0.,6) 

COLL NUMEER C 220 . ,260. ,-4. ,FM0D,0. , 1, ' C II ) ' ) 
COLL ZNEWF 

SI RIPETE PER UN OLTRO MODO 
5 CONTINUE 

SI RIPETE PER UN OLTRO GIORNO 
7 CONTINUE 

DI50TTIU0ZI0NE DEL PLOTTER 
COLL ZDEND 
COLL ZEND 

FINE REGOLORE DELL' ONOLISI 
WRITE! 1 ,*) 

WRITE! 1 , *) ' **FI NE REGULORE DELL '' ONOLISI**' 



REWIND ITWPE 
STOP 

L 

C FINE DEL MR IN 
END 
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FTN7X , 5 
SFTLESCO, 1 ) 

C 

PROGRAM FFREGRESS! ) , < 880323 . 1 4 1 5 > 

C 

C * CRERZIONE DI DISC FILES DEI DOT I STRNDORD ***'" ..* 


C * QUEbTO PRQGRPMMR OPERR SLII DRTI STPNDRRD 
C * 

C * 

C * ESSO CRER DISC FILES R PAR TIRE DRI DRTI I MMRGRZZ I NOT I SO NOSTRO 
C * CIObCUN FILE CONTIENE DRTI RELOT I 01 RD UN SINGOLO PR5SRGGI0 
L * RVENDU L I OE ' FISbPTG LE STRZICJNI RTTIVE, IL MODO ED IL GIORNO - 

* DELL'INIZIO DEL PRS5PGGI0 

* IL FILE CONTIENE TUTTI I DRTI NELLE ERNDE 5 E X E LR FRFQUENZfl 

* DELL' 05CIL.LRT0RE . 


* N.B. 51 E' RSSUNiO CHE FRO UN PRS5MGGI0 E L'RLTRO Cl 5IRN0 8 ORE * 

* * 

** X p X X X X * X X X X X * - X X - V- V y v V V v v y y y y V V V ^ y v y y v s/ y y y y y y y y y y y y y ^ y 


DICHIRRRZIONI 

INTEGER*^ LMRT (160,15), I STORE (3, 730) , RNGDRT C 120 , 6 ) 
INTE6ER*4 I YERR , I YD, ZERO 
DIMENSION ID5RC6) , IDSTC6) , I MODE ( 4 ) 

CHRRRCTER*3 NY 

COMMON IPO, IPOS, I POX, IRUTO, I STORE , RNGDRT , LU 


ESEGUIBILE 

I NFORMRZ I ON I DI I NGRE5S0 / USC I TO 
WRITEd , 100) 

100 FORMRT (IX, 20 / , 80 ' * ' , / ) 

WRITE (1 , *) ' CRERZIONE DI DISC FILES DEI DRTI STRNDORD' 

WRITEC1 !*) ~ w * 

WRITEd ,*) 

WRITEd , *) ’ QUE5T0 PROGRRMMR CRER DISC FILES fl PRRTIRE DRI' 
WRITEd,*) ' NR5TRI STRNDORD.' 

WRITEd,*) ' CIRSQJN FILE CONTIENE DRTI RELRTIUT RD UN SINGOLO' 
WR-l 1 1( 1 ,*) ' PHSSRbGIO NELLE DLIE BRNDE S E X E LR FREOUENZR' 
WRITEd,*) ' DELL' ' 05CILLRT0RE . ' 

WRITEd ,200) 

200 FORMRT (1 X , / , 80 ' * ' , / ) 

10 WRITEd, *) ' **vrPRRRMFTRI PER L'* RNRLI5I : ' 

WRITEd ,'( "NOSTRO IN LETTURR (803)#: ")') 

RERD( 1 , *) I TOPE 

I F ( I TOPE . NE . 8 . PND . I TOPE . NE . 9 ) GO TO 10 
REWIND I TOPE 

WRITE! 1 , ' ( "QURNTI RECORDS? : " ) ' ) 

RERD( 1 , *) NFI 

WRITEd ,'( "OUOI L ' ' RNPLISI RUTOMRTICR? <1=SI>: " ) ' ) 



RE0DC1,*) IOUTO 

WRITEC 1 , ' C " bTOMPONTF tt : ")') 

READC1,*) l-U 
IFCLLJ.NE.'I) THEN 
WR I TEC LU ,300) 

300 F0RMATC5X, ’FILE NAMES' , 6X , ' IN DAY FIN_DAY , 6X , 

+ ' YEAR IN MIN FIN MI N ' , 6X , ' MODE TRb Dbb RFC DbS', 

+5X , ' R GbL ’ , 7X , ' DEC ’ , 6X , ' RTLT ' ) 

WRITECLU,*) 

END IF 

20 WRITEC 1 , ' C "GIORNCl INIZIHLE? < DA 1 M 3b5>: ")') 

READM,*) IFIOYI 

IF ( I DRY I . LT . 1 . OR . I DRY I . GT . 365 ) bO TU 20 
30 WRITE Cl , ' C "DELL' 'ANNO: '13 ")') 

READ Cl,*) I YEAR 

I F C I YEOR . FT . 3S . OR . I YEAR . LT . 1 ) 50 TO 30 
ZERCJ = 0 
NR = 0 
NT = 0 
NM = 0 

DO 1 1=1,6 
NR=NR+1 

WRITE Cl,*) ' **** 0 PER TERMINARE ' 

WRITEC 1 ,' C "5TAZI0NE RICEVENTE? <12,14,42,43,61,63): ")') 

READC 1 , *) IDSRCNR) 

IFC IDSRCNR) .EQ.O) THEN 
NR=NR-1 
60 TO 40 

IFC IDSRCNR) .NE. 12. FIND. IDSRCNR ). NF . 14 . RND . IDSRCNR) .NE.42.flND. 
+ IDSRCNR) . NE . 43 . AND . IDSRCNR) . NE . 61 . AND . IDSRCNR ) . NE . 63) THEN 
NR=NR- 1 
WRITEC 1 ,*) 

WRITEC1,*) ' ****5TAZI0NE NON RICONOSCIUTA !***** ' 

END IF 

1 CONTINUE 
40 DO 2 1=1,4 

nm=nm+i 

WRITE Cl,*) '**** 0 PER TERM I FARE' 

WRITEC 1 , ' C "MODO? < Dfl 1 0 4>: ")') 

READ Cl,*) IMODECNM) 

IFC IMODECNM) .EQ.U) THEN 
NM=NM- 1 
GO TO 50 
END IF 

I F C I MODE C NM ) . LT . 1 . OR . I MODE C NM ) . 6T . 4 ) THEN 
NM=NM- 1 

WR I TE Cl,*) ^ 

WRITEC 1 ,*) '****MODO NON RICONOSLlUTQ ! *****' 

ENDIF 

2 CONTINUE 
SO DO 1=1,6 

NT=NT+ 1 



n n nn 


WHITE Cl,-*) '**** 0 PER TERMINARE' 

WRITFC 1 , ’ C “STOZIGNF TRA5MI TTENTE? < 12 , 1 4 , 4? , 43 , 6 1 , 63 > : 11 ) ' ) 
READ Cl,*) IDSTCNT) 

IFC IDSTCNT) .EQ.O) THEN 
NT = NT - 1 
60 TO 51 
ENDIF 

IFC IDSTCNT) .NE. 12. AND. IDSTCNT) .NE.14.flND. IDSTCNT) . NE . 42 . RND . 
+ I DST C NT ) . NE . 43 . RND . I DST C NT ) . NE . 6 1 . RND . I D5T C NT ) . ME . 83 ) THEN 
NT = NT- 1 
WRITE Cl ,*) 

WRITECI,*) ’ ****STRZIONE NON RICGNOSCIUTfl ' 

ENDIF 
END DO 

LEGGE I DRT I RNGOLRRI 
51 WRITECI,*) 

WRITECI,*) ' LETTURR FILE DAT /FILE#1 ' 

OPEN C 70, FILE= ‘ DAT /FILE#'! J , I05TRT = 1050 , STATUS = ‘ OLD ' , ERR = BO) 
RERDC70 , *) RNGDRT 
CLOSE C 70 ) 

GO TO 70 

60 WRITECI,*) 'ERRQRE #',1050,' IN OPEN IL FILE DRT /FILE#'! ' 
WRITECI, 'C" VUOI ABORT I RE IL PROGRAMMA? <Y/N>: ")') 

READ Cl, 400) NY 
400 FORMAT CR3) 

IFC NY . EQ . ' Y ' . OR.NY.EQ. 'YES' .OR. NY.EQ. ' 1 ' ) 

+5T0P ’ ***F I NE IRREGOLARE DELL'’ ANALI5I***' 

GO TO 10 

51 FI55RN0 I PARRMETRI 

70 DO 5 1 1 = 1 , NR 
ID55R = ID5RC ID 
ID55T = IDS5R 
13=0 

DO 5 12=1 ,NM 

71 I MODES = I MODE C 1 2 ) 

I F C I MODES . GT . 2 ) THEN 

13=13+ 1 

IFC. 13. GT. NT) GO TO 5 
I D55T = I DST CIS) 

IFCID55T.EQ. ID55R) GO TO 71 
ENDIF 

WRITECI ,*) 

WRITECI,*) 'RICEREA: MODO =',IM0DE5,' D55 RICEVENTE= ' , IDS5R , 
+' D55 TRA5MITTENTE= ' , IDS5T 
E 

E I N I Z I AL I ZZAZ I ONE 

REWIND I TAPE 
IP0 = 0 
IP0S=0 
IPOX=G 
I BTF=0 



I3Tfl=G 
i n i t = n 
NEWDRY =0 
I DRY = I DRY I 
DO 3 11 = 1,3 
no 3 JJ= 1,730 
I STORE C 1 1 , J.J)=0 
3 CONTINUE 

I Y D = I YEAR* '1 000 . + I DRY 
ITC=480 


C 

C LE 6 GF 


1 000 


UN RECORD 
WKITEC 1 ,*) 

WRITEC 1 ,*) ' LEGSENOu DRL NflSTRU' , I i Mrfc 


NRL =0 _ , 

CRLL XTRPEC'I , I TREE , LMRT , 4800 , LEN) 

NRC = NRC+ 1 

I F ( LEN . EO . - 1 ) GO TO 2000 

IF CL EN.NF. 48 Q 0 ) WRITE Cl,*) ' **ERRORE 


IN 


C PER CIP 5 QJNR RI 6 R DEL BUFFER 
90 CONTINUE 

DO 4 I = 1 , 1 bO 


C RICERCE DEL 6 I 0 RN 0 RILHIE 5 T 0 
IGO=LNRTC I , 1 ) -I YD 
IFCIGO.LT.O) bO TO 4 
IFCI 60 .GT. 0 ) THEN 
I DRY -I DRY + 1 
I YD= I YD + 1 

IF ( I NIT . EQ . 0) bO TO yO 
NEWDRY =1 
XBTR=1440- IBTF 
I STORE ( 2 , 9 ) = IBTF 
I BTF=0 
END IF 


C 

C RNRLISI DEI DRTI 

IST=LMflTC I , 4) / 10000 
IRM=LMRT C I , 4 ) - 1 ST * 10000 
ISR= IRM/ '100 
I M = I RM - 1 SR* 1 00 -10 
IH = LMRT C 1 , 2) /1 0000 
I MS = LMRT C 1 , 2 ) - 1 H * 1 0000 
IBT= IH* 60 + IMS / 100 


C 

C FINE 


DEL PRSSRGblO 

IFC C IBT-IBTF*IBTR) . bT . I TO 
CRLL CRER 
NEWDRY =0 
I BTF= I BT 
IBTR =0 
ENDIF 


THEN 


L-ETTURfi : 


l.UN= ' , LEN , '**’ 



c 

C DOTO DO TR05FERIRE NEL BUFFER DI L.ISC I TO 

IFC I5R.EQ. ID55R.0ND. IM.EQ. IM0DE5.0ND. IST.EQ. IDSST) THEN 
IPG=IPO + 'l 
C 

C PRINO DOTO PER II. BUFFER DI U5CITO CRIEMPIMENTO DELLO TE5T0TQ) 
IFCIPO.EQ.1) THEN 
I F C I OUTO . NE . 1 ) THEN 
WR I TE ( 1 , * ) 

WRITEC 1 , *) ' DOTI PER GUESTO INPUT' 

END IF 

I STORE C 1 , 8 ) = I DOY 
I STORE C 1 , 3 ) = I DOY . 

I STORE (2,8)=IBT 
I STORE C3, 1 ) =2 
I STORE C 3 , 2 ) = I YEAR 
I STORE C 3 , 3 ) = IDSST 
I STORE C 3 , 4 ) = I D55R 
I STORE (3,5) = INODES 
I N I T = 1 
NEWD0Y=0 
END IF 
C 

C IMMOGOZZI NOMENTO DEI DOTI 
I BTF= I BT 

I STORE Cl , S ) = I DOY 
I STORE C 2 , 9) = IBT 
I T = IBT- I STORE C2 , 8) +1 1 
IFC NEWDOY . EQ . 1 ) IT=IT+1440 
IFCIT.GT.721) THEN 

WRITE Cl , ' ("ULIOI OBORTIRE IL PROGROMMO? <Y/N> : " ) ' ) 
REODC 1,400) NY 

IFCNY.EQ. ' Y' . OR . NY . EQ . ‘ YES' . OR . NY . EQ . ' 1 ' ) 

+STOP '** FINE IRREGOLORE DELL ' ' QNGLISI***' 

ENDIF 

I STORE C 1 , IT) =LMOT C I ,6) 

IST0REC2, IT) =LNOT (1,7) 

I5T0REC3 , IT) =LMGT C I ,3) 

C 

C GGG I ORNQMENT0 DEI DOTI NELLO TESTOTO 
I F ( LMOT C I , 6 ) . NE . ZERO ) THEN 
IFOS= IP05+ I 
I STORE C 1 , 2 ) = J DOY 
I STORE C 1 , 4 ) = IBT 
I STORE C 1 , 7) = IT 
IFCIP0S.EQ.1) THEN 
I STORE C 1 , 1 ) = IDOY 
I STORE C 1 , 3) =IBT 
I ST ORE ( 1 , 6) = IT 
ENDIF 
ENDIF 

IFC LMOT C I , 7) . NE . ZERO) THEN 
IP0X=IP0X+1 



I STORE C 2 , 2 ) = IDQY 
I STORE C 2 , 4 ) = IBT 
I STORE C 2 , 7) = IT 
IFCIPOX.EQ.D THEN 
I STORE ( 2 , 1 ) = I DAY 
I STORE (2 ,3) =IBT 
1ST ORE ! 2 , S) = I T 
END IF 
END IF 
END IF 
C 

C 51 R I PETE PER UN'OLTRU RIbO 
4 CONTINUE 

C 51 RIFETE PER UN OLTRQ RECORD 

IFCNRC.NE.NFJ) bO TO 1000 


C FINE RICERCO DEI DOT I 
2000 COLL CREO 

I YE OR'! = LMOT ( 1B0 , 1 ) / 1000 . 

IDPY1 =LMOT C'l 60, 1 ) - 1 YEAR 1*10C0 . 

WRITE C 1 , *) 

WRITE!!,*) 'FINE OL GIGRNU: ' , ID0Y1 , 
WR I. TEC 1 , * ) 


C 51 PROCEDE CON OLTRI POROMETRI 
IF C 13. 61. 0) 60 TO 71 
5 CONTINUE 


DELL' ' ONNQ ' , I YEAR 


L 

c 


OPZIONE PER RIPETERE L’ RNQLISI CON OLTRI DOTI 
IECI0UT0.EQ.1) WRITE! 1 , *) ' ' 

WRITECLU , *) 

REWIND ITRPE 

WRITE! 1 , ' C " VUOI IMM060ZZIN0RE OLTRI DOTI? 
READ Cl, 400) NY 

IFCNY.EQ. ' Y’ .OR.NY.EQ. ' YES' .OR. NY. ED. ' 1 ' ) 


< Y/N> : " ) ' ) 

60 TO 10 


C 

C FINE RE60L0RE DELL'ONOLlSl 
REWIND ITRPE 
WRITE!'!,*) 

WR I TE C "I , * ) ' **F I NE REGOLOKE 
STOP 


DELL' ' ONQLISI** ' 


C 

c FINE DEL MO I N 
END 


L 

c xxxxk&^x-a-axw*********'* 

C QUESTO SUBROUTINE CREO I DISK FILLS 
SUBROUTINE CREO 


INTEGER*^ I STORE! 3, 730) , QN6DQTC 120 , 6) 
CH0R0CTER*14 NFILE 
CH0R0CTER*3 NY 



r ' n n n n n n n 


COMMON IPO, IPOS, I POX, IPUTO , I STORE , PNGDPT , LU 
C 

C 5E Cl SONO DPT I 

IFCIFO.NE.O) THEN 

IF C IPUTO . NE . 1 ) WRITEC 1 , *) 'Cl SONO', IPO,' DPTI ' 

C 

C RTCERCP DEI DPTI PNGOLRRI 
CRLL 5FRRCH 

CRERZIONE DEL FIL.ENPME 

WRITEC Nr I LE , 1 00 ) I STORE (3,2) , I STORE (3,4), I STORE (1,6) 

100 FORMAT ( ' DAT /ROD' , 12,12, 13) 

SALTO DELLE RICHIESTE DR INPUT SE IL LPUORO E' AUTOMATICO 
IFC IPUTO. EQ. 1) THEN 
NY = ' YES' 

GO TO 10 
END IF 

SCRIVE IL DI5C FILE RICHIE5T0 
WRITEC 1 , *) 

WRITEC 1,*) 'IL FILE ' , NFILE , ' CONTIENE I SEGMENT I DPTI:' 

IFC I STORE (3, 1 ) . EQ. 2) WRITE Cl ,*) '***VRLORI REGRESS' 

WRITEC 1 , *) ' ***BPNDE : S, X E FREQUENZR' 

WRITE ( 1 , *) ' ***GIORNO I NIZIRLE : ' , I STORE (1,8) 

WRITEC 1 ,*) '***GIORNO FINALE ISTOREC 1 , 9) 

WRITEC1 ,*) '***PNNO: ' , I STORE (3, 2) 

WRITEC1,*) '***MINUTO I NIZIRLE: ' , I STORE (2, 8) 

WRITEC 1 , ■*) ' ***MINUTQ FINALE: ' , I STORE (2, 3) 

WRITEC 1 ,4) '***MODO: ' , IST0RE(3,S) 

WRITEC 1,*) ' ***DSS TRASMITTENTE : ' , IST0REC3 , 3 ) 

WRITEC 1,*) '***DSS RICEUENTE : ' , I5T0REC3, 4) 

WRITEC 1,*) '***RSCENSI0NE RETTR : ' , I STORE C 3, 7) 

WRITEC 1 ,*) '*#*DECLI NAZI ONE: ' , IST0REC3,8) 

WRITE Cl ,*) '***RTLT: ' , I STORE C 3, 6) 

WRITE C I ,*) 

WRITEC 1 , ' ( " OK RD I MMRGPZZ I NRRE? <Y/N>: ")') 

RERDC 1,200) NY 
200 FORMAT C A3) 

IFCNY.EO. 'A' ) STOP '***FINE IRREGOLARE DELL' ' RNRLI5I***' 

SCRITTLJRR DEL. FILE SE RICHIESTQ 
10 IFCNY.EO. 'Y' .OR. NY. EQ. 'YES' .OR. NY. EQ . ' 1 ' ) THEN 

IFCLU.NE.1) WR I TE C LU , 300 ) NFILE , ISTOREC 1 , S) , ISTORE C 1 , 3) , 

+ ISTORE C 3 , 2 ) , ISTOREC 2, 8) , I STORE C2, 9) , ISTOREC 3, 5) , ISTOREC 3, 3) , 
+ I STORE C 3 , 4 ) , I STORE (3,7), I5T0REC3,8) , I5T0REC3,S) 

300 FORMRT C 1 X , PI 6 , 1 1 C I 1 0 ) ) 

WRITEC 1 ,*) 

WRITECI,^) 'SCRITTURR DEL FILE ', NFILE 
ISTOREC 1 ,5)=IP05 
ISTOREC 2 , 5) =IPOX 
ISTORE C3 , 9) = IPO 

OPEN C 70 , FILE^NFILE , I OSTRT - 1050 , FRR = 20 ) 



WRITE (70,*) I STORE 

CLOSE (70,1 05TRT = 1 05C , ERR =30 ) 

ELSE 

WRITE! '1 , * ) 'FILE ' , IMF TLE , ’ NONCRERTO' 
ENDIF 
C 

C RFINIZIRLIZZRZIONF DEI BLlFFEEb 
IP0 = 0 
IP05=0 
IPOX=0 

no i n = i,3 
no 1 JJ= 1 , 730 

I ST ORF ( 1 1 , • U ) = 0 
1 CONTINUE 
ENDIF 
C 

C RITORNO RL MR IN 
RETURN 


C 

C ERRORR IN OPEN FILE 70 

20 WRITE (I,*) '**EKRURE IN OPEN FILE 
STOP 


, NFILE , ' #' , I ObO , '**' 


L 

C ERRORE IN CLOSE FILE 70 _ __ kicrTl _ , 

30 WRITEC1,*) ' **ERRORE IN LL.ObE FILE , NFILE, tt 

STOP 


i osl , '** 


c 

c 

c 


FINE Dl CRER 
END 

GUEST A SUBROUTINE CERCR I 
SUBROUTINE SEARCH 


DRTI RN60LRRI 


INTEGER-x-4 I STORE ( 3 , 730 ) , RNGDRT ( 1 20 , 6 ) 

INTEGER*4 5UMMIN , MINMED , IDHM , ~DHM 
CHRRRCTER*3 NY _ 

COMMON IPO , IPOS , I POX , I AUTO , ISTORE , RNbDR l , LU 

r 

C R I CERCR DEI DRTI RNGOLRRI 

SUMMI N= ISTORE (2,8)+ ISTORE ( 2 , S ) + 1 . 

I F ( I STORE C 1 ,8) . NE . I STORE ( 1 ,3) ) blJMM I N = 5L1MM I N + 1 440 . 

M I NMED = SUMM I N / 2 

IF (MINMED . 6T . 1440) M I NMED =M I NMED -1440 
D4= 10000 

I DHM= I STORE (1 , 8 ) *D4 + MI NMED / 1 60*1 00 + MOD ( MINMED , 1 60 ) 

KK = 0 

DO 1 K=1 , 120 _ tn 

IF ( I5T0RE(3 , 2 ) . LT . RNGDRT (K , 1 ) / 1000) bO TO 10 
FDHM= ( RNGDRT (K , 1 ) - 1 STORE ( 3 , 2 . ) * 1 000 ) *D4 + PNGDRT (K , 2 ) / 1 00 
IFCFDHM.LT. IDHM) GO TO 1 
KK=K 



GO TO in 
1 CONTINUE 

c 

C DOT I NON TROVATI 

10 IFCKK.EQ.O) THEN 

WRITE Cl,#) ’**DATI ANGOLRRI NON TROVAT I ! ** ' 

IFCLU.NE.l) WRITE C LU , * ) 'DOT I ANGOLRRI NON TRQVATI ' 
IFCIAUT0.NE.1) THEN 

WRITEC I , ' C "QUO I ABORT I RE IL PROGRRMMA? <Y/N>: ")') 
READC 1,200) NY 
200 FORMAT ( AO ) 

IFCNY.EQ. ' Y' .OR.NY.EQ. 'YES' . OR . NY . ED . ' 1 ' .) 

+ biOP ' *** I RREbOl .ARE FINE DFl_ PROGRRMMA***' 

END IF 

I STORE C 3 , 7 ) =ANGDAT C 120,5) 

I STORE C 3 , 6 ) = ANGDflT C 120,4) 

I STORE (3,9) = ANGDAT C 1 20 , 3 ) 

RETURN 
END IF 

k:kk=kk-i 

IFCKKK.EQ.G) THEN 
I STORE C 3 , 7 ) =RNGDRT C 1 ,5) 

I STORE (3,8) = ANGDAT C 1 ,4) 

I STORE (3,9) = ANGDAT ( 1 ,3) 

RETURN 
END IF 
C 

C 51 RIEMPIE IL BUFFER 

FS 1 = ANGDAT C KKK , 2 ) - ( ANGDAT ( KKK ,2 5/1003*100 

FH 1 = ANGDAT CKKK , 2 )/ D4 

FM1 = C ANGDAT C KKK , 2 ) -FH 1 *D4 ) / 1 00 

F52= ANGDAT (KK , 2) - (ANGDAT (KK , 2) / 100) *100 

FH2= ANGDAT (KK , 2) /D4 

FM2 = ( ANGDAT ( KK , 2 ) -FH2*D4 ) / 1 00 

ITS =30 

ITD= IDHM/D4 

ITH= ( IDHM-ITD*D4) / 100 

ITM=CIDHM-ITD*D4-ITH*100) 

IFCFH2-FH1 .LT.O) THEN 
FH2=FH2+24 

IF ( I DHM-FDHM .LT.O) ITH=ITH+24 
ENDIF 

T2=F52/60. +FM2+FH2*60. 

T 1 =F5 1 /SO . +FM1 +FH1*60 . 

T= I TS / 60 . + I TM+ I TH*60 . 

DT = T2-T 1 

IFCDT.EO.O) THEN 

WRITEC!,*) ' **R I PET I Z I ONE NEI DATI ANGOLRRI**' 
IFCLU.NE.l) WRITECLU,*) ' RIPETIZIONE NEI DATI ANGOLRRI ' 
IFC IAUTO.NE. 1 ) THEN 

WRITEC 1 ,'( "VUO I ABORT I RE IL PROGRAMMA? <Y/N): ")') 

RE ADC 1,200) NY 

IFCNY.EQ. ' Y' .OR.NY.EQ. 'YES' .OR.NY.EQ. ' 1 ' ) 



n n 


RNRL I bl *-*-*' 


+ 5T0P ' **FINE IRREGOLPRE DELL'' 

ENDIF 

I STORE C 3 , 7 ) =RNbDRT ( KK , 5 ) 

I STORE C 3 , 8 ) =RNGDPT C KK , 4 ) 

I STORE (3,3) = RNGDRT (KK , 3 ) 

ENDIF 

S i =RNGijRT i. K K i 5 ) -RNGDRT (K.KK , S J 
I STORE ( 3 , 7 ) =ST / DT* (T-T'l ) + RNGDRT ( KKK , S ) 

ST = RNGDRT (KK, 4 ) -RNGDRT (KKK , 4 ) 

I STORE ( 3 , 8 ) = ST / OT-* ( T - T I ) + RNGDRT ( KKK , 4 ) 

ST = RNGDRT (KK , 3) -RNGDRT (KKK , 3 ) 

I STORE ( 3 , 6 ) = ( 5T / DT* ( T - T 1 ) * RNGDRT ( KKK , 3 ) ) / 1 0 
C 

L RITORNO RLL_R CHIMMRTR 
RETURN 

FINE DI SEARCH 
END 



TAB . 1 


REGRESS MEASUREMENT STRUCTURE 


1 TIME 1: UTC (seconds 1950) 

2 TIME 2: composite word containing UTC coded as year, day - 

of year, hour, minute, second and fractions of second-; 
these information are contained in the decimal part of time 2 

as : 

YYDDDHHMMSSf f f f f 

3 IDENTIFICATION WORD: composite word containing sampling rate 

in hundredths of seconds in the first seven decimal figures; 
the 8th figure represents the downlink band <1=5, 2=X> , the 
9th one the tracking network (DSN-1), the 10th and 11th the 
transmitting station, the 12th and 13th the receiving 
station, the 14th the data type (l=Doppler) , the 15th the 

tracking link mode (2=two way, 3=three way), the last- 
figures are not used; the scheme is: 

1 . SSSSSSSBNTTRRKW 

4 OBSERVABLE: observed frequency shift in Hertz, (same quantity 

as in the ATDF tapes) . 

5 FREQUENCY: oscillator frequency in Hz. The transmitted 

frequency is obtained multiplying FREQUENCY by 96. 

6 PASS IDENTIFIER: composite word containing spacecraft 

i dent i f i cat i on code . 

7 COMPUTED: frequency shift predicted by ODP 

8 RTLT : round trip light time in seconds 

9 HOUR ANGLE: composite word containing the spacecraft hour 
angle in degrees, minutes, seconds and fractions of seconds 
in the first figures; the code is: 

1 . DDMMSSff f f 

10 DECLINATION: composite word containing spacecraft declination 
in degrees, minutes, seconds and fractions of seconds; the 
code is the same as in item 9. 

11 AZIMUTH: composite word containing spacecraft azimuth in 

in degrees, minutes, seconds and fractions of seconds; code 
as in item 9. 



12 


ELEVATION AT RECEIVING TIME: 
s p a c e c r a f t eleva 't i o n a t l- h e 
minutes , seconds and fractions 


composite word containing 
receiving time in degrees 

of seconds; code as in item 9; 


13 El E'VATION AT TRANSMITTING TIME : composite word containing 

spacecraft elevation at the transmitting time; code as in 

item 9, 


14 DERIVATIVE 0F r THE ELEVATION AT RECEIVING TIME: in deg . /s 

15 DERIVATIVE OF THE ELEVATION AI TRANSMITTING TIME: in deg . /s 

16 RESIDUAL: Doppler residual in Hz, computed from OBSERVABLE 

minus COMPUTED (item 4 minus item 7) . 


17 REJECTION CODE: if this code is 0 the data are ok. 

18 WEIGHT: data quatity indicator (used by ODF ) . 

19 CALIBRATION: media calibration to Doppler residuals (in Hz). 

20 CALIBRATED RESIDUAL: calibrated residual in Hz. It is computed 

by subtracting: RESIDUAL minus CALIBRATION (items 16 and 19). 



STRUCTURE OF RECRESS RECORDS 


rows -> 

1 

6 

9 

33 

57 

o 1 1 j m n s I 1 

• 

N1 

H2 

TIME TAG 1 

TIME TAG 2 

ID. WORD 

i o 

Hi 

H2 

OBSERVABLE 

FREQUENCY 

PASS ID. 

V 

3 

HI 

H2 

COMPUTED 

RTLT 

HOUR ANGLE 

4 

HI 

H2 

DECLINATION 

AZIMUTH 

ELEV . T3 

If) 

HI 

N2 

ELEV. T 1 

ELEU. DER, 

T 1 ELEC. DER. 

6 

HI 

H2 

RESIDUAL 

REJECTION 

WEIGHT 

~7 

/ 

HI 

H2 

CAL IE?.. COMP. 

MODELS 

TIME TAG 1 


TAB . 3 


EXAMPLE OF REGRESS RECORD 


The data start from time 20h 00' 29“ of the DOY 308 year 1980 


1 1 

1 2 - 

1 3- 
1 4- 
1 5 
1 6 - 
1 7 
1 8 
1 9 
1 10 
1 11 
1 12- 
1 13 
1 14 
1 1 5 


, 9732816299999999994-009 
. 26909466283331 9999 + 00 5 
. 26909022 59070722 5 5+ 00 5 
. 6 1004 72 584 1368 102!)- 001 
, 544776425522538426+002 
. 443692624774428168+000 
, 000000000000000000 
. 803082000299999998+016 
, 220277400000000000 + 008 
. 101894300547065213 + 00 5 
. 23635580744061 1173 + 003 
, 28382 4 0 5 4 044811 2 3 9 - 0 0 2 
, 000000000000000000 
. 0 0 0 0 0 0 0 0000 0 0 0 0000 
, 100060001114141200+017- 


, 803082000299999998 + 016 
. 220277400000000000 + 008 
. 101894300547065213 + 005 
.236255807440611173+003 
. 2836240540448 1 1 23 9 “002 
, OOP r o o o 0 OOP n n n n n n n 
, 000000000000000000 
. 1 00C 60002 1 14121 300 + 0 1 7- 
, 103080310204800000+017“ 
. 4 10241 0374896060 19+002- 
, OOOGGOOOOOOOOOOOOO 
, 16494637177069948 5-003“ 
. 0000 000 00000000G00 
, 973281 629999999999+009 
. 269156408999980000 + 005 


. 100060001114121300 + 017 
. 103080310204800000+017 
. 410241037489606019 + 002 
, 000000000000000000 
, 164946371770699485-003 
000000000000080000 
. 973281629999999999+009 
. 986680276999990000 + 005 
. 986664161639321474 + 00 5 
. 61004 7258413681014-001 
. 544776425522538426 + 002 
. 161 15360668 5252271 + 001 
, OOOOOGOOOOOOOOOGOO 
. 803082000299999998 + 016 
, 220277400000000000 + 008 


o m i s s i s 


1 94 
1 95 
1 96 
1 97 
1 98 
1 99 
1100 


, 000000000000000000 
, 100060001114121300 + 017- 
.1030803102048 0 0 0 0 0 + 017“ 
. 4 15250900514604227+002- 
. 0 0 ( j 0 (3000 0 0 00000 0 0 0 
. 1 13775830962588772-003- 
. 000000000000000000 


973281749999999999+009 
268720277666650000+005 
2687 15821903272889+005 
609566732206611566-001 
5 4 4 9 4 3 5 9 22872944 70 + 002- 
445576337711003134+000 
000000000000000000 


803082002299999998+016 

220277400000000000+008 

101894286495735869+005 

236819925387191006+003 

285344413867476923-002 

000000000000000000 

000000000000000000 



STRUCTURE OF STANDARD RECORDS 


1 


4 

5 

6 

7 

8 
9 

10 
11 
1 2 
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YYDDD 

HHMMSS 

FREQCY 

TTRRNM 

ELupELdw 


year and day 

hour, minute and second 

oscillator frequency in Hz 

trans . and rec . stations, network and mode 
elevation in up and down link in degrees 


S-RES 

X-RES 

S-3/1 1*X 

S-OBS 

X-OBS 

REJCOD 


S-barid residual in mmHz 
X-band residual in mmHz 
quantity used in plasma calibration 
S-band observable (in Hz) 

X ~ b a n d o b s ervable ( i n H z ) 
reject ion code 


WEIGHT datum weight 
not used 


14 


not used 


15 


not used 



TAB . 5 


DATA TRANSFER FROM REGRESS FILE TO STANDARD FILE 


INPUT DATA 


OUTPUT DATA 


1 TIME 1 

2 TIME 2 


3 ID. WORD 


4 OBSERVABLE 

5 FREQUENCY 

6 PASS ID. 

7 COMPUTED 

8 RTLT 

9 HOUR ANGLE 

10 DECLINATION 

11 AZIMUTH 

12i 

[elevations 

13- 1 


not transferred 
r 1) YYDDD 

2) HI-IMMSS 

t- band indicator 
L 4 ) TTRRMM 

9) S-OBS or 10) X-OBS 

3) FREQCY 

not transferred 
n o t t r a n sferre d 
not transferred 
n o t t r a n s f e r red 
not transferred 
not transferred 


5) ELupELdw 


[EL. DERIVATIVES 
1 5 * 

16 RESIDUAL 

17 REJECTION CODE 

18 WEIGHT 

19 CALIBRATION 

20 CAL. RESIDUAL 


not transferred 

6) 5-RES or 7) X-RES 

11) RE J COD 

12) WEIGHT 

not transferred, 
not transferred. 



TAB. 6 


EXAMPLE OF STANDARD RECORD 


Data start from time 20h 00' 29" of the DOY 308 year 1980 

1 80308 200029 22027740 141213 54473792 -443692 -1611536 -4182 

-26909 -98668 00000 

2 80308 200029 22027740 141412 54473792 -447566 -1624393 -4549 

-26915 -98690 00000 

3 80308 200029 22027740 144213 54472920 -450149 -1638823 -3197 

-31147 -114208 0 0000 

4 80308 200029 22027740 144313 54472920 -446804 -1639150 236 

-31147 -114208 00000 

5 80308 200129 22027740 141213 54483776 -442122 -1610828 -2805 

-26890 -98599 0000 0 

6 80308 200129 22027740 141412 54483776 -446126 -1623059 -3473 

-26896 -98621 00000 

7 80308 200129 22027740 144213 54482936 -447568 -1638656 -661 

-31130 -114146 00000 

8 80308 200229 22027740 141213 54493752 -445576 -1602485 -8534 

-26872 -98530 00000 

9 80308 200229 22027740 141412 54493760 -447826 -1620098 -5981 

-26878 -98553 00000 

10 80308 200229 22027740 144213 54492960 -446743 -1636628 -389 

-31114 -114085 00000 

. . . omissis 

160 80308 204029 22027740 141412 53653088 -445258 -1612198 -5567 

-26197 -96058 00000 


1 



EXAMPLE OF REGRESS TABLE 


LTi fn 3 ~ ^ as N C 3 CO ? F J « sorirt 

if. 'j\ ro pj f J .'J m f ^ -h <h -h <r ~h — i 


O 

CD 

03 

vi 

C_ 

(0 

03 

> 

d) 

r 

-p 





QC 

C 

LJ 


a. X CL x Q- CL 

CO O C O a 

GO C Q Q O 


CL X CL X CL X 

Co Coda 

c* a a coo 


X X 

C O 
Q 


X X CL X 

Coda 
o a a c 


TA&. 7 



ANGULAR DATA STRUCTURE 


TIME TAG 
TIME TAG 
RTLTxlO 

R Ah 100 
DECK 100 
ELOhIOO 


(high part) , coded as YYDDD 
(low part), coded as HHMM5S 

it is given in seconds with an approximation 
of ,1 sec; the stored value is IFIX (RTLT*10 . > 

the right ascension is given in degrees with 
an approximation of .01 degrees; the stored 
value is IFIX (RA* 100 . ) 

the declination is given in degrees with an 
approximation of .01 degrees; the stored va- 
lue is IFIX (DEC-xlOO. > 

the elongation is given in degrees with an 
approximation of .01 degrees; the stored va- 
lue is IFIX (ELOx-100 . ) 



TAB. 9 


LOGICAL FORMAT OF DOPPLER DATA FILES (ROD) 


rows > 


c o 1 1 j m n s ! 1 

1 

i 

i n 
i c. 

i 

i 

V 3 S 
4 


5 tt 


1 

S INITIAL DAY 
S FINAL DAY 
INITIAL MINUTE 
S FINAL MINUTE 
of 5-band DATA 


X INITIAL DAY 
X FINAL DAY 
X INITIAL MINUTE 
X FINAL MINUTE 
tt of X-band DATA 


3 

DATA ID. ( 1 =A 2=R) 
YEAR 

TRANSMITTING DSS 
RECEIVING DSS 
MODE 


6 LOCATION OF THE 
FIRST S-BAND ITEM 

7 LOCATION OF THE 
LAST S-BAND ITEM 

8 INITIAL DAY 

9 FINAL DAY 

10 0 

11 ! 

S-BAND DATA ! 

, mmHz ! 


LOCATION OF THE 
FIRST X-BAND ITEM 
LOCATION OF THE 
LAST X-BAND ITEM 
INITIAL MINUTE 

FINAL MINUTE 

0 


X-BAND DATA 
mmHz 


RTLT 

RIGHT ASCENSION 
DECLINATION 
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1) S and X stand for S-band and X-band respectively 

2) The value of the item DATA IDENTIFIER is 1 for ATDF data and 
2 for REGRESS data 



TAB. 10 


FILE DAT/ROD 14 14309 


1 

309 

309 

2 

2 

309 

309 

80 

3 

903 

903 

14 

4 

1347 

1347 

14 
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384 

384 

2 

6 

1 1 

11 

10187 

7 

455 

455 

18571 

a 

309 

903 

-22 

9 

309 

1347 

384 

10 

0 
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0 

i l 

-1716015139 

2147483647 

22027744 

12 

2147483647 

2147483647 

22027744 

13 

2147483647 

2147483647 

22027744 

14 

-458095 

-1670888 

22027720 

15 

-457467 

-1671956 

22027720 

16 

-457156 

-1669117 

22027720 

17 

-453513 

-1655412 

22027720 

730 

0 

0 

0 



FIGURE CAPTIOHS 


Fig. 1 The Deep Space Network stations. 

Fig. 2a Flow chart of the program DECODEREGRES. 

2b Zoom of block A of the previous flow chart. 

Fig. 3 Flow chart of the program STATREGRES. 

Fig. 4 Output of the program PLOTREGRES. 

Fig. 5 Flow chart of the program PLOTREGRES. 


Fig. 6 Flow chart of the program FFREGRES. 
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